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Abstract 

We study the instanton effects of the ABJM partition function using the Fermi gas for- 
malism. We compute the exact values of the partition function at the Chern-Simons lev- 
els k = 1, 2, 3, 4, 6 up to iV = 44, 20, 18, 16, 14 respectively, and extract non-perturbative 
corrections from these exact results. Fitting the resulting non-perturbative corrections 
by their expected forms from the Fermi gas, we determine unknown parameters in them. 
After separating the oscillating behavior of the grand potential, which originates in the 
periodicity of the grand partition function, and the worldsheet instanton contribution, 
which is computed from the topological string theory, we succeed in proposing an analyt- 
ical expression for the leading D2-instanton correction. Just as the perturbative result, 
the instanton corrections to the partition function are expressed in terms of the Airy 
function. 
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1 Introduction 



Partition functions are the most fundamental quantities in both quantum systems and sta- 
tistical systems. Although it typically serves only the role of normalization in the physical 
interpretation, it encodes important thermodynamic information such as the free energy. In 
field theories, the partition function, in general, contains the information of not only pertur- 
bation but also instanton effects as well as some global analytical structures. Therefore it is 
natural to study it as a first step. 

M-theory is supposed to unify all our understandings of non-perturbative effects in the 
string theory. In the seminal paper [1], the low-energy effective theory on iV multiple M2- 
branes on the geometry C 4 /Z fc was proposed to be the Af = 6 supersymmetric Chern-Simons 
theory with gauge group U(N) k x U(N)^ k and bifundamental matters Ai :2 , B li2 forming the 
superpotential W = \{A X B X A 2 B 2 - A X B 2 A 2 B X ). 

It was found in [2, 3, 4] that the infinite-dimensional path integral of the partition function 
on S 3 and the expectation values of BPS Wilson loops are reduced to a finite dimensional 
matrix model by the localization techniques [5]. 

There is much progress in this matrix model including the derivation of the degree of free- 
dom iV 3 / 2 [6] expected from the AdS/CFT correspondence [7] and its perturbative completion 
by the Airy function [8]. Recently, the ABJM partition function was rewritten as that of a 
Fermi gas system [9, 10, 11], where the standard method in statistical mechanics is applicable 
and the above properties are reproduced easily. 

The Airy function and the Fermi gas formalism imply a deep structure of M-theory. The 
Airy function has the integral representation of the exponentiated cubic term, which looks 
similar to the partition function of the Chern-Simons theory, and appears also in other contexts 
of M-theory [12]. Therefore we would like to study the partition function further and to extract 
some physical information out of it. 

In our previous work [13], we have proposed a method to compute the exact partition 
function Z k (N) at k — 1 and indeed computed it up to N = 9. Here with a generalization 
of our method and a technique from [14], we perform the computation of the exact partition 
function at k = 1, 2, 3, 4, 6 up to JV = 44, 20, 18, 16, 14 respectively. 

After the computation of the exact values, we proceed to study the non-perturbative effects. 
We can extract the non-perturbative corrections from the exact values of the partition function 
by subtracting the perturbative result. The main motivation of this paper is to explore an 
analytic form of the non-perturbative correction. It is known that there are two kinds of 
instantons that induce the non-perturbative corrections to the ABJM partition function. One 
is called the worldsheet instanton, and the other the D2-instanton (membrane instanton). 
From the viewpoint of the gravity dual, the Type IIA string theory on AdS 4 x CP 3 , the 
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worldsheet instanton comes from the fundamental string wrapping the holomorphic cycle 
CP 1 C CP 3 [6, 15], while the D2- instanton comes from the D2-brane wrapping the Lagrangian 
submanifold MP 3 C CP 3 [16]. The worldsheet instanton behaves as e ~^V^ while the D2- 
instanton as e - w ^ 2kN [ n the large N limit. 

Let us summarize our results here. Using the Fermi gas formalism (and some results of 
the topological string), we can know the expected forms of the instanton corrections to the 
partition function. In the Fermi gas formalism, it is useful to consider the grand partition 
function. The grand potential (the logarithm of the grand partition function), in general, 
receives the non-perturbative corrections as 

- w - 4 pert V) = 4 D V) + jr^) + • • • , (i.i) 

where \x is the chemical potential, and jj^ CTt \fi), J® 2 (/j,) and J^ s (fj,) are the perturbative, D2- 
instanton and worldsheet instanton contributions, respectively. We note that the dots in (1.1) 
represent an additional oscillatory contribution and, if any, bound states of the worldsheet 
instanton and D2-instanton. The oscillatory behavior will be studied in subsection 3.2 in 
detail, which we ignore here. The D2-instanton and the worldsheet instanton corrections are 
expected to take the following forms: 



■TO = E J f (m V) = E [4 m V + & i m V + 4" 

m=l m=l 

oo oo 

71=1 n=l 

We would like to determine the prefactors of those instanton corrections. The prefactors of 
the worldsheet instanton corrections are determined by the results of the topological string on 
local F , where the first few prefactors are 

dW = - d {2) = - - d {3) = - + - (13) 

k sin 2 (f)' dk 2sin 2 (f) sin 2 (f)' 3 sin 2 (f ) + sin 2 (f ) " ^ 

On the other hand, there are no systematic ways to determine the prefactors of the D2- 
instanton corrections. 1 In this paper, we propose the analytic form of the leading D2- instanton 
correction, 

(i) 4 f nk \ l(i) 2 (ilk 

ai = 7— cos — , oi — —7— cos — 

k n 2 k V 2 J k 7rtan(f ) V 2 



c (1) 



2 5k k 1 



3A; ' 12 ' 2sin 2 (f) 7rtan(f, J 



cos ( ^] . (1.4) 



1 One approach is to consider the 'WKB' (small k) expansions of the prefactors as in [9]. However, it is 
hard to determine the higher order corrections in this way. 
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We arrive at these expressions from many constraints (by matching the exact data, the nu- 
merical data from TBA, pole structure, small k behavior etc.) and a few assumptions. We 
believe that these results are valid for any k. We also find the higher instanton corrections at 
k = 1,2, 3, 4, 6 from the numerical fitting. Here is the summary of our results: 



(np) 



V + n + 1/4 



71^ 



52/x 2 + ///2 + 9/16 



2tt 2 



+ 2 



-8/i 



+ 
+ 



32 



,-12/x 



736/^ 2 - 152///3 + 77/18 
3^ 

2701/x 2 - 13949/i/48 + 11291/768 



+ 466 



Ji np) = 



4/i 2 + 2/i + 1 



71 z 



7T 



e-2M + 



52/i 2 + /i + 9/4 



T-(np) 
J 3 



(np) 



J, 



(np) 



+ 
+ 

4 

-e 3 
3 

e^ + 
4 

-e 3 



736/i 2 - 304/i/3 + 154/9 
3^ 



2tt 2 



32 



+ 2 



e -16 M + e ,( e -20^ 



-6m 



270V - 13949/i/24 + 11291/192 

4/^2 + ^+ 1/4 20" 
3^ + ¥ 



+ 466 



-8m 



+ 0(e 



-10|U\ 



2e~^ + 
4/x 2 + 2/x + 1 



-4/x _ ^ 

9 



27T 2 



2e~3^ + 



e -2, + 16 e -3, + 
O 



52/i 2 + /x + 9/4 



47T 2 



+ 2 



e' 4fl + 0(e~ bfl ), 



\ 



4/i 2 + 2/x + 1 20 
3^ + ~9~ 



9 



(1.5) 



Once the instanton corrections to the grand potential are fixed, we can easily translate 
them into the instanton corrections to the partition function, 



7 (exact) 



(AO 



4 pert) (iV) 



1, 



(1.6) 



where Z k peTt \N) is the perturbative result in [8], whose explicit form is given by (3.18). The 
leading corrections of both instantons are written by the Airy function (and its derivative): 



Z^(N) = (C?(N + 2 - B k )a^ + <£>) ^ 
v J Ai 



lA Ai[C fc 1/3 (iV + 2-ff fc )] 



[C k 1/3 (N-B k )} 



_ c _ 1/%{1) M[C k l/ \N + 2-B k )] 
k k M[C~ 1/3 (N-B k )] 
z ws (N) = d(1) Ai[C^\N + l-B k )] 
k M[C- 1/3 (N-B k )] ' 



(1.7) 
(1.8) 
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with 

Some remarks on our results are mentioned in order. 

Firstly, we emphasize that our results include all the 1/JV corrections. One can check that 
in the large N limit, both instanton corrections show the correct exponentially suppressed 
behavior as mentioned above. 

Secondly, the worldsheet instanton correction (1.8) is divergent at k — 1,2, and the D2- 
instanton correction (1.7) is also divergent for even k? Therefore we cannot apply the above 
result to these values of k. In these cases, the D2-instanton and worldsheet instanton correc- 
tions are mixed, and the sum of them totally become finite. At k — 1, 2, for example, such 
mixed instanton corrections are given by 

Z™+™(N) - (2N + C ' + 29 ) A ^i 1/3 (^+f)] <#* Ai'[Cr 1/3 (* + f)] (1 10) 

Zl {n) -K 2N+ y + t) Ai[ cr 1/3 (iv-§)] i"Ai[cr 1/3 (^-|)]' ( } 



nwW/ x / \ A[[C 2 1/3 (N + 1)} 2/3 Ai'[Gr 1/3 (7V + |)] 

Z D2+ws , N s = Un + C 2 + 7) 1 2 - TT r — - 2C 2 2/ — -^ttt — ■ (1-H) 

Ui[C- l/3 (N-l)} Ai[C" 1/3 (iV-i)] 



One can check that the leading large N behavior of (1.10) precisely reproduces the result in 
[14]. We note again that our results include all the 1/N corrections. 

Finally, for k > 2, the leading non-perturbative contribution comes from the worldsheet 
instanton. Since our result is valid for any iV and k, we can take the 't Hooft limit. Our 
worldsheet instanton correction (1.8) must give the correct answer in this limit. We check 
that the genus- zero and genus-one worldsheet instanton corrections to the free energy studied 
in [6] are indeed reproduced from our result (1.8) after taking the 't Hooft limit. This is a 
non-trivial check of our result. 

The paper is organized as follows. In section 2, we review the method to compute the 
partition function exactly. We generalize our previous method in [13] to arbitrary level k. 
In section 3, using the exact partition function obtained by this method, we numerically 
study the behavior of the instanton corrections. We focus on the grand potential and the 
partition function. We find that the grand potential shows an oscillatory behavior and that 
it is explained by the periodicity of the grand partition function. By the numerical fitting, 
we determine the instanton corrections to the grand potential at k — 1, 2, 3, 4, 6. In section 4, 
we compute the worldsheet instanton corrections by using the results in the topological string 
on local F . As a non-trivial test, we take the 't Hooft limit of our results, and confirm 



2 For odd k, the D2-instanton correction (1.7) trivially vanishes, and the leading correction starts from 
0(e" 27rv ^). 
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that our results correctly reproduce the worldsheet instanton corrections to the genus-zero 
and genus-one free energies up to 3-instanton. In section 5, we propose an analytic form of 
the D2-instanton correction, and check its validity. Section 6 is devoted to conclusions and 
discussions. In appendix A, we prove a new relation found in this paper. In appendix B, we 
summarize the exact values of the partition function at k — 1, 2, 3, 4, 6. 



2 Exact Results of Partition Function 

In [9] the ABJM partition function 



Z k (N) 



V!) 2 / 
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(iV!)2 J (2n) N (2n) N 

is rewritten in terms of the partition function of a Fermi gas system 

d N q 

o-e5jv 

with the density matrix given by 

11 1 1 
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(2.3) 



After introducing the chemical potential /i or the fugacity z = e M , the grand partition function 

oo 

3*M = 1 + E ^(^e^, (2-4) 

AT=1 

is given by 

{-z) n 



Sfc(/i) = det(l + zp) = exp 



E 

n=l 



n 



■Trp r ' 



(2.5) 



This means that if we want to know the partition functions for various N, all we have to do 
is to compute Tr p n for n = 1, • • • , N. 

As explained previously in [13], it is difficult to compute Trp n because of the complexity 
in convoluting the density matrix. This difficulty was, however, overcome by a lemma from 
[17] stating that when a kernel K(qi, q 2 ) takes the form 



K(qi,q 2 ) 



E{qi)E{q 2 ) 
M( qi ) + M(q 2 y 



(2.6) 
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we can compute its power K n (q 1: q 2 ) by 
where (p e (q) is defined by 

0%) = / ^K\q,q>)E(q>) . (2.8) 

We can compute (p e (q) recursively 

d>\q) = ^jf <i)E{q')¥~\q), (2-9) 

with the initial condition <f)°(q) = 1. Though we can apply this method to the density matrix p 
directly, there is a more convenient way to compute the grand partition function as explained 
in [13]. 

Let us briefly review the method found in [13]. As our density matrix is invariant under 
the parity 

-92) = p(qi, <h), (2.10) 

we can decompose it into two parts: 

I \ l \* I \ l \ - P(9i>92) ±p(qi,-q 2 ) , 011 v 



Since p + (and also p_) is written as the form (2.6) 

£(gi)£(g 2 
cosh(y) + cosh(^- 

with 



P+(gi,g2) = ^xT^Tn I ^u^v ( 2 - 12 ) 



. , cosh(^) 

= m Xm > (2 - 13) 

^/2/c cosh(|) 



we can compute p\ from a series of functions <f> (q) as 



n-l 



= co 3h ( ¥ )?I-^co S h( f ) E(-l)V(^-'te). (2.4) 

The important fact is that using (2.9) the functions £ (<?) are determined by the following 
recursive relation: 

1 f°° dq> 1 cosh(f^) 



tffq) = / — ; ( 2 .15) 

VW 2A;cosh(^)y_ 00 27r2cosh(^)cosh(^) V ^ 1 K ' 
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One can show that the grand partition function is factorized as 



H(/x) = det(l + zp + ) det(l + zpJ). 

In [13], we found that the density matrices enjoy an interesting property. 

det(l + ^-; ) y , , 
det(l-zp + ) 

Thus the grand partition function is rewritten as 

oo 

S(/i) = det(l-zV + )-^^(0)/, 



(2.16) 



(2.17) 



(2.18) 



£=0 



and determined only by the even parity functions (f> (q). 

In [13], we further found that the eigenvalue problem of the density matrix p is mapped 
to that of an infinite-dimensional parity-preserving Hankel matrix, whose form is given by 
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Thus the grand partition function is also expressed as 

E(fi) = det(l + zp) = det(l + zH). 



(2.20) 



In this paper, we find a general formula for arbitrary parity-preserving Hankel matrices, which 
reduces to (2.17) for our current case. Let H± be the even/odd parity block of the original 
Hankel matrix, 
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Then it satisfies the following relation 

det(l + 



H 



(hi 


h 2 


h 


\ 


h 2 


h 


hi 




hs 


h± 


h 5 




\ 






■) 



(2.21) 



det(l - zHj 



\ e o| 



zH 4 



\v), 



(2.22) 



7 



with 



|e ) = (l,0,0,---) T , b> = (1, 1, 1, - - -) T - (2-23) 

In appendix A, we shall give a proof of (2.22) for a certain class of Hankel matrices. 

Using these methods we computed the partition functions at k — 1,2,3,4,6. The exact 
values are given in appendix B. As can be imagined from [18], the partition functions have 
close relations with tan(|). Therefore it is especially clean for these values of k. 

3 Numerical Study 

As summarized in appendix B, we obtain the exact values of the partition functions at k = 
1,2,3,4,6 up to N max = 44,20,18,16,14, respectively, following the procedure in section 2. 
After computing the exact partition function, let us extract important information from it. 
Here we numerically study the non-perturbative corrections to the partition function. We shall 
plot various physical quantities computed by our exact values and fit them by the theoretically 
expected forms [9, 6, 16]. 

3.1 Grand Potential 

First, we start with the grand potential Jfc(/x) defined by 

JM = !og . (3.1) 

From the theoretical expectation the prefactor of the instanton corrections are at most quadratic 
functions of fi (see (1.2)). In fact, we find numerically that the grand potential Jfc(//) (k = 
1, 2, 3, 4, 6) generally can be expanded as 3 

oo 

JM = 4 PCTt) M + E , 4 n) = 4"V + /t V + 7i n) • (3.2) 

n=l 

In the following we shall plot the grand potential and their instanton corrections, fit them by 
the quadratic functions and read off the coefficients numerically. 

However, not all of the coefficients are read off from this analysis of Jk{^)- As we shall 
explain below we will encounter some oscillations in fj, and it is difficult to study the exact 
values of the coefficients due to large errors caused by the oscillations. Some of the values are 
actually first read off from the analysis of the partition function Z k (N) as in subsection 3.3 
which does not have the oscillating behavior, while some of them are found after we understand 

3 We expect that this property holds for fceZ. However, we cannot establish it at present. 
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Figure 1: (a) The perturbative, (b) 1-instanton, (c) 2-instanton and (d) 3-instanton corrections 
of the grand potential Jk{fJ>) at k — 1. The dots represent the grand potential obtained by the 
exact partition function while the solid lines represent the fitted functions. 

the oscillating behavior well as in subsection 3.2 and subtract its contribution. It is simpler, 
however, to present our result as if we read all of the exact values from the grand potential 
and explain various behaviors later. Hence, we shall take this path in presenting our results. 

Let us start to study the grand potential at k = 1. In figure 1(a), we have plotted 
the grand potential Ji(/i) (for selected values of /i) in dots together with its perturbative 
expression J-[ per ^(/z) in a solid line. The grand potential Jk(l^) (3.1) is given by the grand 
partition function Sfc(/x) (2.18) with the summation of N truncated at our maximal values 
A^max, while the perturbative expression is given by 

Here the constants Bk and Ck are given by (1.9), and Ak was conjectured in [9, 19] as 

Ak . log ^ + 2a _ X) 1 T d 1 (3 _ 1 _ 3 N 

»7r 6 k 3 J e KX — 1 \x 6 x xsinh x/ 

We find a perfect match, especially for large /i. Therefore, we can plot (//) = (Ji(/i) — 
j[ pcTt \fi)) /e~ 4/i next in figure 1(b) and fit the result by a quadratic function. We find the 



behavior matches well with 



We continue this process further and plot j[ 2 \y) = (J(fj) — j[ pert \y) — e ^J^\y))/e 8m in 
figure 1(c). Now the plot starts to oscillate. We can fit the result by the sum of a quadratic 



function with the oscillations cos(^) and sin(^), 

= - 5V+ <t'\ +9/16 + 2 + jr c (^ ), 

7f»- W = 2cos(^-|). (3.6) 

Similarly, the 3-instanton jf\y) = (■Ma*) - ^^V) ~ e~ 4/1 Ji\y) ~ ^ jf\y))/e^ is 
plotted in figure 1(d), and the fitted function results in 

J? = 736^-15^3 + 77^18 _ 32 + ^ 

, f >« ( ^-i ( 8, + 1)si „(^-!|)-3 2 oo S (^-I|). (3.7) 

Note that our plot range of /i changes slightly. This is because in going to higher and 
higher instanton effects, we have more and more numerical errors and it is difficult to find a 
sensible plot for large y. 

We can repeat the same analysis at k — 2,3,4,6 since we have exact values. The results 
computed from the exact values we find are plotted in figures 2, 3, 4, 5, respectively together 
with the fitted functions. The results are summarized in (1.5). Using the results in (1.5), 
we can study various instanton effects. This will be done in the subsequent sections. Before 
proceeding to it, let us give a physical interpretation for the oscillating behaviors and present 
the numerical study of Z(N) without the oscillations. 

3.2 Oscillations in </&(//) 

In this subsection, we would like to explain the oscillation in Jfc(/i) from the periodicity of 
the grand partition function. Since the grand partition function H fc (/i) = e Jfc< ^ is a sum of 
e Nfl defined by (2.4), it must be a periodic function of \x with the periodicity 2ni. However, 
the naive expression of derived from the method of statistical mechanics in small k 

perturbation is not invariant under the shift + 27r«. 

To remedy the invariance, we define the 27ri-periodic grand partition function by 



= J2 e J £ aivc ^ +2 ™). (3.8) 
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Figure 2: (a) The perturbative, (b) 1-instanton, (c) 2-instanton and (d) 3-instanton corrections 
of the grand potential Jfc(/x) at k — 2. The dots represent the grand potential obtained by the 
exact partition function while the solid lines represent the fitted functions. 



Here J£ aive (/i) means the grand potential obtained from the naive small k, large [i asymptotic 
expansion, namely J£ aive (/i) = J^ p (//) + (//) with jj^ crt \fi) given by (3.3). Through 
this periodic extension (3.8), the grand potential acquires the oscillatory part 



jm = jr ivc (/i) + 4 osc (a0 

where J£ sc (/i) is given by 

TO = iog(i + -^£e-<r'» 



(fj,+2nin) 



(3.9) 



(3.10) 



For instance, the oscillatory term in the worldsheet 2-instanton comes from the ±2ni shift 



of the perturbative part jjf CTt \fi). Since 



e 4 pcrt V) 



exp 



8/i 
"T 



we find the oscillatory term 



I' 

e k 



(pert) 



(P) 



2 cos 



3k 



2n[C k fi 2 + B k - — 



e k 



(3.11) 



(3.12) 
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Figure 3: (a) The perturbative, (b) 1-instanton, (c) 2-instanton and (d) 3-instanton corrections 
of the grand potential J k {fj) at k — 3. The dots represent the grand potential obtained by the 
exact partition function while the solid lines represent the fitted functions. 



in the worldsheet 2-instanton correction. Note that this exactly reproduces (3.6) at k — 1. 

We can proceed further to study the oscillatory term in the worldsheet 3-instanton for 
k > 2 from the worldsheet 1-instanton contribution 



•CM 



7(P ert )/ . „ ■ \ I ,WS(1), . „ .s ,(pcrt), „ .n ,WS(1), n -\ 
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e k + 4d k sin — sin 

K 



/ 14 
2n[C k ^ + B k - — 



12m 

e *> . 



(3.13) 



For = 1,2, however, since the D2 m-instanton j^ 2(m ) ^ e" 2m ' 1 (m = 2, 1 for = 1,2) has 
the same contribution as the worldsheet 1-instanton term {ij) ~ e £ , we have to be 

careful. The mixing between these two kinds of instanton effects will be studied in more details 
in section 5. Using and in (3.2), we find that the oscillatory term in the worldsheet 
3-instanton for k = 1, 2 is given by 



J ( 3 )osc (/i) = _ Alx (2a k l) ii + sin 2n(c kl x 2 + B, 



3A; 
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Figure 4: (a) The perturbative, (b) 1-instanton, (c) 2-instanton and (d) 3-instanton corrections 
of the grand potential Jfc(/i) at k = 4. The dots represent the grand potential obtained by the 
exact partition function while the solid lines represent the fitted functions. 



?7r 2 'a^ cos 



(3.14) 



Note that this expression for k = 1 reproduces (3.7). 

As we have seen in the previous subsection, we actually observed an oscillatory behavior 
in the numerical analysis of Jfc(/i). Since the amplitude and the phase of oscillation perfectly 
matches (3.13) (and (3.14)), we strongly believe that our numerically observed oscillation is 
originated from the 27rz-periodic extension of the naive grand partition function. 

One important consequence of the periodicity of the grand partition function is that in 
the integral transformation from Sfc(/i) to Zk(N) we can extend the integration range from 
[—711,711} to (— zoo,ioo) 



Zk (N) = / ^.e^M-"" = / -^V* , M-"a* . (3.15) 
J -iri 2tt2 J -too 2ttz 

Since the last expression of (3.15) does not involve J£ sc , there is no oscillation in the canonical 

partition function Zk(N). The last expression in (3.15) also explains the good agreement of 

the Airy function with the exact partition function Zk(N). 4 



l To obtain the Airy function expression, we need to deform the integration contour from [—iri, iri] to 
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Figure 5: (a) The perturbative, (b) 1-instanton, (c) 2-instanton and (d) 3-instanton corrections 
of the grand potential </&(//) at k — 6. The dots represent the grand potential obtained by the 
exact partition function while the solid lines represent the fitted functions. 

3.3 Partition Function 

In the previous subsections, we have studied the numerical behavior of the grand potential 
Jfc(ju) and observed an interesting oscillating behavior. Though interesting from a theoretical 
viewpoint, the oscillating behavior causes a larger error in finding the exact values of instanton 
prefactors. As mentioned above, this behavior will disappear after an integral transformation 
from the grand potential to the canonical partition function. For this reason, as a complemen- 
tary analysis, let us look at the partition function as well. We find exactly the same values 
from the fitting of Jfe(/i) and the fitting of Zj.(N). This agreement gives a strong support for 
the validity of our method. 

Let us consider the case k — 1. In this case, the non-perturbative correction to the grand 
potential starts from 0(e~ 4fl ) as found in [13]. The (naive) grand potential is generically given 
by (3.2). We would like to determine the constants a[ , and 7^ from the exact data. 

(— ioo,ioo). Naively, this deformation causes the error 0(e~ 2M */ fc ) = 0(e~ 7r V 2Ar / fc ), where /i» is the value of 
the saddle-point, as was explained in [9]. However, such a correction does not actually appear in the difference 
of the exact partition function and the perturbative result. The procedure here naturally explains this fact. 
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Since the partition function is computed as 



Zl (N) = f °° ^Vr ivc w-^ 

J-ioc 27rt 



-L 



100 

too 



[l + (4V + /3i> + 7i ij )^ + 



2tt? 



(3.16) 



the 1-instanton correction normalized by the perturbative partition function is given by 



Z ( ^{N) 



1 



z[ pert \N) J-i 



^ e ^ rt V)-^ (Q; ( V + ^(i) + 7 (D )e -4, 



z (pert) (iv) 



2vri 



(cr 1 (iV + 4- J B 1 )aS 1) +7! 1) )- 



lA Ai[Cr 1/3 (7V + 4-.B 1 )] 

Aiicr 173 ^-^)] 



_ r -i/3^(i) Ai / [gr 1/3 (iV + 4-i? 1 )] 
1 * Ai[C^\N- Bl) ] ' 

where Z^ pcrt ^(A^) is the perturbative contribution, which is given by [8, 9], 

4 pcrt) (iV) = C~ 1/3 e A * Ai \C- 1/3 (N - B k 



(3.17) 



(3.18) 



The expression (3.17) is the form expected from the Fermi gas consideration. We compare this 
result with the non-perturbative correction computed from the exact values of the partition 
function using our definition of z[ np \N) in (1.6). To fix the constants (5^ and 7^, we 
fit the exact data Z^ p) (N) / for 20 < N < 44 by the expected form z{ 1] {N)/e~ 2irym , 
and find that af \ f3^ and are very close to (in about 14-digit agreement) 



a 



(i) 



7T Z 



ft 



(1) 



1 



7T Z 



7i 



(i) 



4tt 2 



(3.19) 



This result, of course, is consistent with (3.5) obtained from the fitting of Ji(^i). In figure 6(a), 
we show the exact data z[ np) (N)/e- 27rVm (the dots) and our result z[ 1] (N)/e~ 27TVm with 
(3.19) (the solid line). Our Airy function expression of the instanton correction (3.17) is in 
excellent agreement with the exact data from A" = 1 to A" = 44. In the large A" limit, Z^\N) 
is asymptotically expanded as 



Z[ 1] (N) = e 



-27TV2JV 



2Ar- 137r2 - 2 y^+ 1697r4 + 1167r2 + 8 
2V2n 32tt 2 

7r(21977r 2 + 4012) 28561tt 4 + 73424vr 2 - 32768 



384V2iV 



+ 



12288N 



+ 0(N- 3/2 ) 



(3.20) 
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(a) (b) 

Figure 6: The 1- and 2-instanton corrections to the partition function at k = 1. We plot (a) 
Z\ iP ^ /e~ 27TV ^ and (b) — z[^) / 'e" 4lV ^ against N. The instanton corrections in terms 

of the Airy function (the solid lines) excellently agree with the exact values (the dots) even 
for small N. 

The leading term 2Ne~ 2nv ^ agrees with the result in [14]. Note that our result is exact at 
all orders in 1/JV. 

Once we fix a[ , pjp and 7^ , we can proceed to determining ot[ \ p± and jx . The 
2-instanton correction to the partition function is given by 

2vr 4 aS 2) C 1 + 3d + 244 



z[ 2 \n) 



2N Z 



-N 



128 



- Cf /3 I N 4 



4tt^{ 2) 




32 



327r 4 7i 2) + l„ 2 122vr 4 af } + 119 3721^ Ai[C^ 1/3 (N + 8 - Bj] 
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Ai^/^JV-SO] 



Ai'[Cr 1/3 (A + 8- Bt)] 



16 ~' ' 8 J Ai[C^ 1/3 (N ~ B 1 )} 
In the same way as above, from the numerical fitting we find 



(3.21) 



a 



(2) 



26 



ft 



(2) 



1 

'4^2' 



7i 



(2) 



9 



+ 2 . 



(3.22) 



7T^ 47T* 327T 2 

These are in about 14-digit agreement with the fitted values. In figure 6(b), we show the exact 
data {Z[ I1P \N) - Z^ ] (N)) / e~ 4nVm and our result z[ 2) {N)/e~ AnVM with (3.22). 

Repeating this method, we find the values of ol^\ P± and 7^ up to n = 4. Furthermore 
we can apply this procedure to other levels. Our results at k — 1, 2, 3, 4, 6 are summarized in 
(1.5). 

Using (1.5), one can easily know the instanton corrections to the partition function. We 
have indeed checked that the instanton corrections to the partition function obtained from (1.5) 
precisely agree with the exact data. The results in (1.5) show that the prefactors at k = 2, 6 
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are strongly related to those at k — 1,3, respectively. Currently, we do not understand the 
physical origin of this relation. It would be interesting to study what this means. 

4 Worldsheet Instantons 

In the previous section, we have obtained the exact non-perturbative correction from the 
numerical study of the grand potential Jk(fj) and the partition function Zk(N). The non- 
perturbative correction consists of several contributions written as (1.1). In this section, we 
would like to determine the prefactors of worldsheet instanton corrections. The prefactors 
of the D2-instanton corrections will be considered in section 5. The prefactor d k n ^ of the 
worldsheet n-instanton is expected to be independent of /i. We can determine d k n ^ from the 
topological string on local F (note that F = P 1 x P 1 ). 

4.1 Topological String on Local Fq 

As discussed in [20, 6], the ABJM matrix model (2.1) (and its generalization to the ABJ model 
[21]) is related to the lens space matrix model [22, 23] 5 by an analytic continuation of the rank 
of the gauge group, and it is further related to the topological string on local F via a large 
iV duality. Using this chain of dualities, we can compute the worldsheet instanton corrections 
in the ABJM theory from the knowledge of the topological string on local F . 

In general, the worldsheet instanton sum in the topological string has the following struc- 
ture [25] 

n=1 g,d 

where \ s denotes the string coupling of topological string theory. 

The free energy of ABJM matrix model and the free energy of topological string on local F 
are related by a Fourier transformation [9] due to the wavefunction nature of the topological 
string partition functions [26]. Since the canonical partition function Z k (N) is obtained from 
the integral transform of grand partition function e Jk ^\ it is natural to identify the grand 
potential Jfc(/i) with the free energy F top (Ti, T 2 ) of topological string on local Fo- Here T\ and 
T 2 denote the Kahler parameters of the two P^s in F . Since the ABJM slice is defined by 
T\ = T 2 , we arrive at the relation 

J k (fi) = F to ^T 1 = T 2 = T), (4.2) 

5 The ABJ(M) matrix model and the lens space matrix model are also closely related to the supermatrix 
model in [24] . 
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with T and \i related by [9] 



(4.3) 



To compare the topological string and ABJM theory, we should be careful about the 
normalization of string coupling. The string coupling A s in the topological string theory and 
the string coupling g s = ^ in the ABJM theory are related by 



where [9] 



Therefore, A s is given by 



n 2sr-2ptop _ \2g-2ptop 



Ac 



47T 

9s = ~~k ' 



and the worldsheet instanton sum (4.1) becomes 



*" op =££»H 2 
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(4.4) 



(4.5) 



(4.6) 



(4.7) 



The Gopakumar- Vafa invariants n 9 ? of local F were computed in [28] . Up to total degree 



d\ + d,2 = 3, the GV invariants n 9 - are given by 



l i,o ~ n o,i ~ 2, n° ! — 4, n^i — n 12 — 6, 
and n 9 -. — (g > 1) for the above cases. Hence F top of local F is expanded as 



(4.8) 



^to P _ _ 



n=l 



e -nTi , e -nT 2 g-n(Ti+T 2 ) 3( e -n(2Ti+T 2 ) , e ~n(T 1 +2T 2 )\ 

H 1 + 



2nsin 2 2p 
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(4.9) 



Using the relations (4.2) and (4.3), we can read off the coefficient of the worldsheet n- 
instanton e~~£~ . The first three coefficients are given by (1.3). More generally, we can write 
down the worldsheet instanton correction to the grand potential in the form 



n,g,d=l 



■2ku\ 29 - 2 (-l)^e-^ 



k J 



n 



(4.10) 



where n g d denotes the weighted sum of GV invariants on local F [27] 



■nr. 
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d\-\-d2= d 
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TV 



di,d,2' 



(4.11) 



For instance, using the data n° A = —48, n\ = 9, n\ — (g > 2) and (4.8), the 4-instanton 
coefficient dfj^ is found to be 

J$ = q _ 1 ]_ / , ! 9 n 

fc sin 2 f 2sin 2 f 4sin 2 ^- 1 J 

In a similar manner, we can determine the higher instanton coefficients using the data of GV 
invariants in [28, 27]. 

4.2 Worldsheet Instanton in the 't Hooft Limit 

In the previous subsection we have determined the worldsheet instanton correction J^ vs (/x) 
from the topological string on local Fo- As a consistency check, we consider the 't Hooft limit 
of the worldsheet instanton corrections derived from the results in the previous subsection and 
compare with the known results of the genus-zero and genus-one free energy. In the 't Hooft 
limit, the free energy has the genus expansion 

oo 

F(g„\) = log Z(g s , A) = ^^" 2 F 9 (A), (4.13) 

9=0 

where g s — is the string coupling and A = y is the 't Hooft coupling. 

One easily finds that the worldsheet instanton corrections to the partition function up to 
3-instanton are given by 6 



Z ™V(N) = d M Ai[C * 1/3{N+ ^- Bk)] 



' M[C- 1/3 (N-B k )} ' 

(4 1} ) 2 \Ai[q: 1/3 (iv + |-5 fc )] 



Ai[C- 1/3 (7V - B k )\ 

lW\3\ A:r/n--l/3/ 



^ ws(2) (iv) = (4 2) + %^) 

Z™®(N) - (d^ + + ^) 3 \ ^ C ^ N + T-Bk)] 

where d^ (n = 1,2,3) are given by (1.3). It is convenient to introduce N = N — |j in the 
analysis below. In the large N limit, 7 Z k NS< " 1 \N) is asymptotically expanded as 



z ws (1)(A0 



f i 5y/2ir rW /25tt 2 1 A k_ a ((ii_\'' 



sin 2 (27r/A;) 



3k 2 V N V 9k4 k2 J N \\N J 



(4.15) 



6 We assume that the worldsheet instanton corrections up to 3-instanton are not mixed to the D2-instanton 
corrections. This assumption is true for k > 6, and the 't Hooft limit, of course, satisfies this condition. 
7 Strictly, we consider the large A = N/k limit (strong coupling limit). 
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If we further expand this in g s = we obtain 



7 WS(1) 



(TV) 



-27TV2A 



-97 A + o 



1 



0{g 2 s ) 



(4.16) 



3 QV2itV^ 4tt 2 A 

where A = y = A — ^. Thus the worldsheet 1-instanton correction to the free energy are 
given by 

~1 5 1 



7 WS(1) 
3=0 



(A) = 



<:? 1J (A) 



3 6^x71 4tt 2 A 



-27TV2A 



(4.17) 



In a similar way we obtain the genus-zero and genus-one part of the 2- and 3-instanton 
corrections from (4.14) as 
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where we have used the relation between the instanton corrections to the partition function 
and to the free energy: 
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Let us compare these results with the result in [6] . One immediately finds that the instanton 
corrections to the genus-zero free energy exactly agree with the result in [6] (see also [19]). 
Let us consider the genus-one free energy. The genus-one free energy is given by 



F g=1 = -log^(r), 



(4.21) 



where t)(t) is the Dedekind eta function, and the modulus r is related to the 't Hooft coupling 
through 
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(a) (b) 
Figure 7: (a) The 2- and (b) 3-instanton corrections to the genus one free energy. 



In the strong coupling limit, the genus-one free energy behaves as 



= ^ - ^log(2V2A) + 



3 QV2ttVX 4tt 2 A 



+ 0(e 



-47TV2A 



). (4.23) 



The 1-instanton correction exactly agrees with (4.17). It becomes messy to compute the 
higher instanton corrections. Instead, we numerically evaluate (4.21), and extract the 2- and 



3-instanton corrections from it. In figure 7(a), we plot {F g= \— F ( 



(pert) 
9=1 



iws(i; 

9=1 



-47TV2A 



against 



A. The dots represent the numerical data obtained from (4.21) while the solid line represents 
our prediction (4.18). Our results show a perfect agreement with the numerical evaluation. 



Similarly, in figure 7(b), we plot (F g= i 
3-instanton prediction (4.19) again agrees with the numerical ones. 



^(pert) 
9=1 



-,WS(1) 
'9=1 



7 WS(2) 
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)/e- 



-67I-V2A 



against A. Our 



5 Conjecture of D2-instanton Corrections 

The coefficients of the worldsheet n-instanton corrections diverge when k G 2n/N. How- 
ever, we confirmed numerically that the grand potential itself is finite at these values of k. 
Furthermore, at k — 2n/m (m = 1,2, . . . ), the worldsheet n-instanton correction is just the 
same order as the D2 m-instanton correction because e _4n/i//fc = e _2mAt . Therefore it is natural 
to expect that D2-instanton corrections are also divergent at these values of k, and the diver- 
gence coming from the worldsheet instanton is canceled by the divergence from D2-instanton. 
The sum of these two contributions leads to a finite result. This story imposes a strong 
constraint on the pole structure of the D2-instanton corrections. 



21 



5.1 D2 1-instanton correction 



Let us consider the e~ 2fl term in jj^ p \/j,). At k = 2n with neZ, the e~ 2M term receives two 
contributions: one from the D2 1-instanton and the other from worldsheet n-instanton. The 
coefficient of worldsheet n-instanton is divergent in the limit k — > 2n. This divergence has the 
general structure 
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with some constant w^ n \ We have checked this behavior up to n 
Vafa invariant of local F in [28] and found 
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Next, let us consider the D2 1-instanton correction 
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The small k expansion of the coefficients a^\b^ and cj^ was computed in 
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From our results at k — 1, 3, we expect that the coefficients of e~ 2 ^ vanish for odd k G Z. 8 Also, 
in order to cancel the poles coming from the worldsheet instanton (5.1), the D2 1-instanton 
coefficients should behave as 
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3 For odd k, it is plausible that there are no D2-branes wrapping 



as in the case of k = 1 discussed 



Further, we have 



,-2irV2kN\ 



in [14]. Thus it is expected that the D2-instanton correction starts from 0( 
numerically checked that this expectation is true at k — 5, 7 by using the TBA-like equations in [13]. 
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From the above conditions, we can fix the forms of af£\b^ and The result is given by 
(1.4). Summing the contributions from the worldsheet n-instanton and the D2 1-instanton, 
the e~ 2fJl term at k = 2n gives the finite answer 



V + 2/1 + 1 



nir 2 



s (n) 



e~ 2 \ (5.6) 



where 



8 W = + ---. (5.7) 
3n 3 y J 



For lower n, we find 



S W = 0, ^=0, s (3) = ^, s (4) = 13, (5.8) 

9 

and (5.6) for n = 1, 2, 3 correctly reproduce our results at k = 2, 4, 6 in (1.5)! 

As a further check of our conjecture (1.4), we compare the instanton corrections to the 
partition function at various k with the numerical results obtained by using the TBA-like 
equations. As explained in [13] (see also [14]), the ABJM partition function can be also 
computed by solving the TBA-type integral equations. 9 In figure 8, for instance, we show 
the instanton corrections at k — |. At k — |, the leading contribution is the D2 1-instanton 
correction. In figure 8(a), the dots represent Zy^ (N) / e~ 7r ^ 2hN computed by the TBA-like 
equations while the solid line represents the D2 1-instanton correction Zf 2 2 (N) / e ~ n ^ 2kN given 
by (1.7) with our conjecture (1.4). In figure 8(b), we plot the next-to-leading contribution. At 
this order, the worldsheet 1-instanton correction is dominant. The dots represent [Z^ 2 \n) — 

Z^ 2 2 {N))/e~ 2lT ^ 2W ^ computed by the TBA-like equations (and (1.7)) while the solid line 
represents the worldsheet 1-instanton correction (N) / e~ 27T V 2N / k given by (1.8). In both 
figures, our theoretical predictions (1.7) and (1.8) perfectly explain the numerical results from 
TBA. We have checked the validity of our conjecture (1.4) for other fc's in the same way. 



5.2 Comment on D2 2-instanton correction 

By the similar reasoning, it seems to be possible to guess the forms of D2 2-instanton coef- 
ficients af\b^ and cf\ However, the constraint is much weaker than the 1-instanton case, 
hence we can not determine the forms of 2-instanton coefficients completely. We need more 
information to fix them. 

9 One advantage to use the TBA-like equations is that we can numerically compute the partition function 
at any k. However, we cannot neglect numerical errors as N grows. 
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Figure 8: In these figures, we plot (a) the leading and (b) the next-to-leading contributions of 
the instanton corrections at k = |. The dots represent the numerical values from the TBA-like 
equations while the solid lines represent our theoretical predictions. The leading correction 
is captured by the D2 1-instanton, and the next-to-leading correction by the worldsheet 1- 
instanton. 



The k — > n limit of worldsheet n-instanton has the general structure 
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We have checked this behavior up to n = 7 and found 
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Naively, these poles should be canceled by the D2 2-instanton. At this order, however, there 
might be another contribution that we will discuss in the next section, i.e., the contribution 
from a kind of bound states of the D2-instanton and the worldsheet instanton. At present, we 
cannot conclude whether such a contribution exists or not. If exists, we need to consider its 
pole structure. We leave this problem as a future work. 

6 Discussions 



In this paper, by generalizing our previous method we continued the exact computation of 
the partition function in the ABJM theory using the Fermi gas formalism. Using these exact 
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values of the partition functions, we further studied the instanton effects. We have found that 
the grand potential contains an interesting oscillatory behavior and it can be explained by 
the periodicity of the grand partition function. The worldsheet instanton effects are under- 
stood from the results in the topological string on local Fo- After subtracting the worldsheet 
instanton, we studied the D2-instanton effects, only little of which has been known so far. We 
proposed an analytic expression of the D2 1-instanton correction which is consistent with all 
of the known properties so far. We would like to discuss some open questions in our analysis. 

In our previous work [13], we have found that the partition functions Z\(N) at k — 1 are 
written as a polynomial of ^ with rational coefficients. Here we have studied the partition 
functions for few other fc's. We find that at k = 2 the partition function are polynomials of 
4j, at k — 3, polynomials of \ and -^=, at k — 4, polynomials of \ and at k — 6, polynomials 
of ^=^. We note in passing a curious observation about the zeros of these polynomials. For 
example, if we replace \ by a variable x in the partition function Z\{N) at k — 1, it becomes 
a polynomial of x of degree [y]. If we order the zeros of this polynomial by their absolute 
values, the first few zeros are very close to and so on. We observed a similar pattern 

of zeros for other fc's. We would like to understand the general structure of Z k (N) better. 

We have determined the analytic expression of the leading D2-instanton correction. We 
would like to study this effect from the worldvolume theory of the Euclidean D2-brane wrap- 
ping on MP 3 . We expect that the prefactor of D2-instanton correction is coming from the 
1-loop determinant of the fluctuation around the classical D2-brane action. To analyze such 
effects, the recent result of reproducing the 1-loop log term from the gravitational computation 
in [29] might be a helpful guide. 

Although the worldsheet instanton coefficients d^ obtained from the topological string 
correctly reproduce the 't Hooft limit of free energy studied in [6], we find some discrepancy 
between the prediction of dfj^ and our results (1.5) obtained by the numerical analysis. For 
example, let us take a closer look at the k — 3,6 cases in (1.5). The coefficients of the first two 
terms in and Jg np ' ) agree with the value of d^ and d^ in (1.3). However, the fourth term 
of and Jq 11 ^ is different from d^ = d^ = —8 predicted from the topological string (4.12). 
Also, for the k — 4 case, the third term of j| np ' ) is different from the prediction d^ = 
while the first term agrees with d^ = 1. We believe that the predictions of worldsheet 
instanton corrections from the topological string (1.3) and (4.12) are correct at generic k, 
and we conjecture that the reason of this discrepancy is due to a kind of bound states of the 
D2-instanton and the worldsheet instanton. This conjecture is based on the observation that 
after a few worldsheet instanton terms which agree with d^\ the discrepancy arises just after 
the first D2- instanton term appeared in . We leave more detailed study of the origin of 
this discrepancy as a future problem. 

Recently, a very interesting relation between the sphere partition functions and the non- 
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perturbative completion of topological string partition functions is proposed in [30]. This 
proposal is reminiscent of the relation between the canonical and the grand canonical partition 
functions (3.15). As discussed in [9], it is natural to identify the grand partition function 
of ABJM theory with the non-perturbative completion of the topological string partition 
function on local F . It would be interesting to study this relation further. 
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A Proof of (2.22) for General Hankel Matrices 

In this appendix, we will prove the identity (2.22) assuming that a Hankel matrix has the 
following integral representation 



We also assume that the weight w(t) is a positive even function of t on the interval (—1,1) 
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(A.l) 



w(-t) 



w(t)>0, fe(-i,i). 



(A.2) 



By the change of variable 



t = t q = tanh -, 



(A.3) 



the matrix element if, 



m,n 



is rewritten as 




(A.4) 



and the trace becomes 




(A.5) 
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where 



2 cosh ■ „ 

n=0 2 



Let us prove the identity (2.22) for this class of matrix 

det(l + zH_) 



det(l - 

where #(2) is defined by 



H(z), (A.7) 



00 1 
If (z) = £(1 - ^ + )oi = (eol^—^lv) . (A.8) 



m=0 



with |e ) and |i>) given in (2.23). 

A.l Proof of the Even Part 

First we consider the even part of (A.7) 

detfl -z 2 H 2 



det(l - z 2 H 2 ) 



H{z)H{-z) . (A.9) 



As discussed in our previous paper [13], using the relation H 2 = H + {\ — \e )(e \)H + the left 
hand side of (A.9) is written as 

det(l -z 2 Hl) . , 1 . . . 

det(l - z 2 H\) = (e °' Y^Hl M ■ 

By applying the lemma of Tracy- Widom to the kernel 

/ /\ PM) +p{q,-q') E{q)E{q') /—— q 

= j = cosh q + cosh q- E(q) = ^(t q )cosh-, (A.ll) 

we find 

21 

The left hand side of (A.12) is rewritten as 

p+ +1 M) = j d ^---j d -^P + {^)---P + M) 

= X/"X! y 2^"' y -^-P2m(q)p2m(qi)---P2n(q2e)p2n(q') 



m=0 n=0 
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= P2m(q)(H 2 + e ) m , n p2n(q') , (A.13) 

m,n=0 

and 4> k (q) in (A. 12) is denned recursively 

<P k (q) = / q'Mq')^- 1 ^') • (A.i4) 



If we expand (p k (q) in terms of P2n{q) 

n=0 wa " 2 n=0 

then al^ satisfies 



fK 00 ,00 



m=0 

The initial condition 0°(g) = 1 corresponds to 

ag> = l (Vm). (A.17) 

In other words, a;^ = ^m- From (A. 15) and (A.16), fc (g) is found to be 

1 00 

0fc(g) = 00^ ^ • (A.18) 

2 ra,m=0 

Setting q = in (A.18), we find 

00 

fe (O) = ^(^) o , m =(e o |^|t;). (A.19) 

m=0 

This implies that if (2;) in (A. 8) can also be written as 

00 

H(z) = 0"(O)2 n . (A.20) 

ra=0 

Finally, taking the limit q, q' — > in (A. 12) we find 

2£ 

«)o,o = ( eo \Hl e \e ) = ^(-l)V fc (O)0 2£ - fc (O) . (A.21) 



k=0 



This is the desired relation (A. 9). 



28 



A. 2 Proof of the Odd Part 

Let us prove the odd part of (A. 8). The kernel p(q, q') can be written as a matrix element in 
the "auxiliary" quantum mechanical system 

P(q,q') = (g|y / ^c pV /^g|g'), [q,p] = 2iri . (A.22) 
In this subsection, we use the notation 

C "=2^h|' S «=2^hf' ^ = tanh |- ( A - 23 ) 
Let us consider the following function 

E^z) = det(l + zR q Tlc q c p ), R q = w^c' 1 . (A.24) 

Using the commutation relation found in our previous paper [13], 

c p c q = -isgllsp, II = 1 - \0 p )(0 q \, (A.25) 

we can see that Si (2) is an even function of z, as follows. We first rewrite Si(z) as 

Ei(z) = det(l + izRgUspU.^ s q ) 
= det(l + izRqSqliSpY^) 

= det(l - zR q c p c q U r ) . (A.26) 
Since Si (2) is a real function for z G R, we find 

E 1 (z) = S 1 (z) t = det(l - z(i? 9 c p c g n t ) t ) = Si(-z) . (A.27) 
On the other hand, S x (z) can also be written as 

Si(z) = det(l + zw(t q )c p ) ■ (0 q \R; \,\., RM 

= det(l + zp) • (O^- 1 -^— ^|0 P ) , (A.28) 

1 + zp 

where 

= \jw{t q ) cosh |, p = sj w{t q )c p sj w{t q ) . (A.29) 
One can easily see that the second factor in the last line of (A.28) is 



Q q \E?-±-E q \%) = 5> B (0)(-*)« = H(-z) . (A.30) 

n=0 

Then Si (—2) = Si(z) implies the following relation 

if(z) det(l + det(l + 



if (-*) det(l - det(l - ' 
This completes the proof of the odd part of (A. 7). 

29 



(A.31) 



B Exact Values of Z k (N) 



Here let us summarize the exact values of Z k (N) (k — 1, 2, 3, 4, 6). Our main task is to solve 
the integral equation (2.15). Here we use the method in [14]. This method allows us to rewrite 
the convolution in (2.15) as the sum of the residues. 



B.l Recursive Relations at Even k 

Let us first consider the even k case. Introducing the variables u = e q ^ k and v = e q '^ k , then 
the integral equation (2.15) becomes 

V? (u) = ~ 7 / dv~, rj-, tzU) (v), (B.l) 

Yy ' 2itu + lJ (u + v)(v k + iy y h y ' 

where f e (u) = <p e (q)- Note that i/ fc ~ 2 )/ 2 is rational when k is even. From the solutions for 
lower £, we observe that the solution <p e (u) have the following form 

£ 

V \u) = J24 ) (u)logiu, (B.2) 

3=0 

where A^\u) is a rational function. We assume that A^\u) has the poles at the roots of 
(u + l)(u k + 1) = 0. We can rewrite the equation (B.l) by using the technique in [14] as 
follows, 

- J__!L_ r dv vk/2 ~ 1 ( v + 1 ) V^MlocH- 
V {U) -2nu + lJ dv ( u + v)(vK + l)^ o A t {V)l ° g V 

1 u J^(2my r v k /*-\v + l) aU) ,. d fhgv 



2tt u + - 

3=0 



1 4i J + 1 A (« + »)(»* + 1) V 2 ™ / 



j=0 J all poles V 7V 7 V 

(-w) fc / 2 (w-l) J^(2iri) j+1 /logu l\ , , N 

= V"^ — 77 — ; — Af (—u)Bj + i I — 1 — ) + fe(u), B.3 

27r(M + l)( M fc + l) ^ j + 1 1 v y J+1 V 27ri 2y ,/n ; ' v ; 

where B n (x) is the Bernoulli polynomial, and fe(u) is a rational function. In the second line 
in (B.3), we have used the integral formula noted in [14], 



Jo 



dv C(v) log' v = j dv C(v)B J+1 , (B.4) 
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where we choose the branch cut of log as the positive real axis, and integration contour 7 is 
chosen as the closed path from +00 to infinitesimally below the cut and then from to +00 
infinitesimally above the cut. Using the identity: 

B n (x + a) = J2 B n - P (a)x p , (B.5) 
p =o 

the first term of the right hand side in (B.3) becomes 

2^ + 1)^ + 1)^1. j + 1 { p ) B ^{2j A ' ( U) {2ni) 

" 2 7 r(r 1 + l)(«* + l)^7TTi j ) B ^-'\2) A > ( U) ) 

(_ u )*/2( u _i) ^ < (2iriy +1 -i fp + 1\ „ /1\ X1 , ,„ , 

jr— p=j-l 
p>0 

We note that B p+1 _,-(l/2) vanishes if p + 1 — j is odd. Thus setting p+l—j = 2m, we obtain 



ri-j + l 

e+i '— 



(» + D(«* + DtJ £r„ ^ + 2 '" I i J^UJ- 4 ' ( " )log *' 



j=0 m=0 

j+2m-l>0 



Comparing the coefficients of log- 5 ' u, we obtain the recursive relations for A^\u). 
For 1 < j < £ + 1, we find 



iH+ii 



(B.7) 

while for j = we find 



4» = - (ti+1)(ufc+1) L — 2^ — 52m U J ^ + ^ (R8) 

v / \ / m=l 

where the rational function /e('u) is given by 

fj u ) = V ^ '- V Res rV-; -^Aj'(v)B j+1 -2- . B.9 

,/n ; 2tt(m + 1)^ j + 1 ^ « ( M + i;)(v fc + l) < v ; 3+1 \ 2ni J V ; 
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where Yl v ^-u means the sum over all the poles except for v = —u. It is convenient to redefine 



iu)f..\ - 1 ( uk/2 v m, 



Ay'(u) = ——- A'\u). (B.10) 

1 v ; (2tt)^'! \u k + I J e K J 



Then we can rewrite the recursive relations as 

r l-j+l i 



_ 1 1 2 i / k/2 \ zm 

m=0 ^ ' 

i ^^2^1 / fe/2 \ 2m 

^>) = -^Em^ ^ 2m - 1) (-«) +/,(«), (B.i2) 

m — 1 ^ ' 



m=l 

where 

_ (-l) m D / I \ .. (-1)"' /' I 

Cm, 



(2m)! B2m U J = 72^)! (,2^ " l ) B ^ (B ' 13) 

u + lj^{ 3 + l)\ v v{u + v)\v k + l) iK>3+ \2m) 

The initial condition is trivially given by 

A^\u) = 1. (B.15) 

The equations (B.11)-(B.15) are the recursive relations to be solved. 

Once these equations are solved, we obtain <p e (u), and thus can compute £ (O) = <^(1) 
and Trp+ n . The matrix element p+ n (wi,w 2 ) is given by 

2n-l 

p^( Ul ,u 2 ) = F k ( Ul ,u 2 ) ^(-VVMvF-^'M, (B.16) 

e=o 

where 

F U, « ^ = («i« 2 )^(«i + l)h + l) m 7 ^ 

fcl Ul ' 2j 2jb( Ul - u 2 )( Ul u 2 - \){u\ + 1)V2( U * + i)i/2- W- 1 ') 

The trace of p+ n is then computed by 

It is easy to see that p+ ra (w, u)/u has the following form 

2n( \ 2n— 1 

*±^=£i&>(„)log'u. (B.19) 
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where E^(u) is a rational function. From the computation for low n's, we observe that (u) 
has the poles at the roots of u 2k — 1 = 0. We assume that this observation is valid for arbitrary 
n. Let u rn be the root of u 2k — 1 = 0, 

/ 7T7TI \ / ITTfl \ 

u m = cos J +isin J , (m = 0,±l,±2,...). (B.20) 
Then, Tr p+ n is computed as 

poo 2™- 1 



2vr Jo 

j=0 J m=0 v 7 

B.2 Recursive Relations at Odd /c 

Next let us consider the odd k case. Changing the variables to t — e q l 2k and s = e q '^ 2k , we 
obtain 

^4^r vry+i/ "' w ' <B - 22) 

where ?/> £ (t) = £ (<?). We observe that the solution takes the form 



^(t) = Y,4 ) (t)lo g n. (B.23) 

3=0 

where a rational function A^\t) has poles at the roots of t 4k — 1 = 0. Let v m be the roots of 
t 4fc - 1 = 0, 

/ 7YVK \ /TTVK\ 

^ = cos(— j+,sin(— j, (m = 0±l,±2,...). (B.24) 
In the same computation in the previous subsection, we find the recursive relations for 



r'' i 



/L,+,(t) - (_1) (f + !)(*»' -D^IJTT I 3 - J S "«-'UJ 

x [4">( It ) + (-l)--Mf(-it)] , (1 < j < [C/2] + 1), (B.25) 



p=0 
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where 

f(t)- 1 * V g^ VRes gfc -V + 1 ) A U) (s)B f^ s \ (B26) 
Mt) -~n¥Ti^ ~7TT~ ^ ££ ( t 2 + s 2 )(s2fc + 1) ^ J • (B.26) 

j=0 m=0 

The matrix element p+ n (ti,t 2 ) is given by 



2n-l 



p 2 -?(h,t 2 ) = F k (h,h) ^(-l)^^)^"- 1 ^^), (B.27) 



£=0 

where 



- (t 1 t 2 ) fc / 2+1 (ti + l)(t2 + 1 ) 

fc(tl ' h) ~ 2k(t* - t 2 2 )(tit 2 2 - \){tf + iy/ 2 (tf + i)V2- ( B - 28 ) 



We further assume that p™(ti,t 2 )/t takes the form 



P?(M) 



n-l 



= E^Wlog^, (B.29) 



3=0 

where a rational function R 2 n{t) has poles at the roots of t 4k — 1 = 0. Thus Tr p 2n is computed 
by 

n-l 



7 />oo " ^ 

«Jo ^ 

j=0 J m=0 v 7 

Let us summarize our procedure here. For even k, we solve the recursive relations (B.ll)- 
(B.15), then compute Trp+ n by (B.21). For odd k, the recursive relations to be solved are 
(B.25)-(B.26). We can compute Tr p+ n by (B.30). These equations do not contain any integrals, 
and we can simply solve by using Mathematica. We have actually solved these equations, and 
have computed Z k (N) up to N = 44, 20, 18, 16, 14 at k — 1, 2, 3, 4, 6 respectively. These values 
are summarized in figures 9-13. 

B.3 A Guess 

The exact values of the partition functions obtained here are all polynomials of 1/ir. At k = 1, 
in particular, the partition function has the following form 

= (B.31) 



=o 
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where the coefficients Z\ (N) are rational numbers. It is, of course, difficult to conjecture 
all the coefficients zf^ (N) for general N. However, we find that the leading and the next-to- 
leading coefficients are simply given by 

Zr(2 M + 1 ,^(-^)% M (^), 

where H n (x) is the n-th Hermite polynomial. One can check that this guess is correct at least 
up to N = 44 by using the values showed in figures 9 and 10. 

References 

[1] O. Aharony, O. Bergman, D. L. Jafferis, J. Maldacena, "N=6 superconformal Chern- 
Simons-matter theories, M2-branes and their gravity duals," JHEP 0810, 091 (2008). 
[arXiv:0806.1218 [hep-th]]. 

[2] A. Kapustin, B. Willett, I. Yaakov, "Exact Results for Wilson Loops in Superconformal 
Chern-Simons Theories with Matter," JHEP 1003, 089 (2010). [arXiv:0909.4559 [hep- 
th]]. 

[3] D. L. Jafferis, "The Exact Superconformal R-Symmetry Extremizes Z," JHEP 1205, 159 
(2012) [arXiv: 1012.3210 [hep-th]]. 

[4] N. Hama, K. Hosomichi and S. Lee, "Notes on SUSY Gauge Theories on Three-Sphere," 
JHEP 1103, 127 (2011) [arXiv:1012.3512 [hep-th]]. 

[5] V. Pestun, "Localization of gauge theory on a four-sphere and supersymmetric Wilson 
loops," [arXiv:0712.2824 [hep-th]]. 

[6] N. Drukker, M. Marino, P. Putrov, "From weak to strong coupling in ABJM theory," 
Commun. Math. Phys. 306, 511-563 (2011). [arXiv: 1007.3837 [hep-th]]. 

[7] I. R. Klebanov and A. A. Tseytlin, "Entropy of near extremal black p-branes," Nucl. 
Phys. B 475, 164 (1996) [hep-th/9604089]. 

[8] H. Fuji, S. Hirano and S. Moriyama, "Summing Up All Genus Free Energy of ABJM 
Matrix Model," JHEP 1108, 001 (2011) [arXiv: 1106.4631 [hep-th]]. 



35 



[9] M. Marino and P. Putrov, "ABJM theory as a Fermi gas," J. Stat. Mech. 1203, P03001 
(2012) [arXiv: 1110.4066 [hep-th]]. 

[10] M. Marino and P. Putrov, "Interacting fermions and N=2 Chern- Simons-matter theories," 
arXiv: 1206.6346 [hep-th]. 

[11] A. Klemm, M. Marino, M. Schiereck and M. Soroush, "ABJM Wilson loops in the Fermi 
gas approach," arXiv: 1207.0611 [hep-th]. 

[12] H. Ooguri, C. Vafa and E. P. Verlinde, "Hartle-Hawking wave-function for flux compact- 
ifications," Lett. Math. Phys. 74, 311 (2005) [hep-th/0502211]. 

[13] Y. Hatsuda, S. Moriyama and K. Okuyama, "Exact Results on the ABJM Fermi Gas," 
arXiv: 1207.4283 [hep-th]. 

[14] P. Putrov and M. Yamazaki, "Exact ABJM Partition Function from TBA," 
arXiv: 1207.5066 [hep-th]. 

[15] A. Cagnazzo, D. Sorokin and L. Wulff, "String instanton in AdS(4) x CP**3," JHEP 
1005, 009 (2010) [arXiv:091 1.5228 [hep-th]]. 

[16] N. Drukker, M. Marino, P. Putrov, "Nonperturbative aspects of ABJM theory," 
[arXiv: 1103.4844 [hep-th]]. 

[17] C. A. Tracy and H. Widom, "Proofs of two conjectures related to the thermodynamic 
Bethe ansatz," Commun. Math. Phys. 179, 667 (1996) [solv-int/9509003]. 

[18] K. Okuyama, "A Note on the Partition Function of ABJM theory on S 3 " Prog. Theor. 
Phys. 127, 229 (2012) [arXiv: 11 10.3555 [hep-th]]. 

[19] M. Hanada, M. Honda, Y. Honma, J. Nishimura, S. Shiba and Y. Yoshida, "Numerical 
studies of the ABJM theory for arbitrary N at arbitrary coupling constant," JHEP 1205, 
121 (2012) [arXiv: 1202.5300 [hep-th]]. 

[20] M. Marino and P. Putrov, "Exact Results in ABJM Theory from Topological Strings," 
JHEP 1006, 011 (2010) [arXiv:0912.3074 [hep-th]]. 

[21] O. Aharony, O. Bergman and D. L. Jafferis, "Fractional M2-branes," JHEP 0811, 043 
(2008) [arXiv:0807.4924 [hep-th]]. 

[22] M. Marino, "Chern-Simons theory, matrix integrals, and perturbative three manifold 
invariants," Commun. Math. Phys. 253, 25 (2004) [hep-th/0207096]. 



36 



[23] M. Aganagic, A. Klemm, M. Marino and C. Vafa, "Matrix model as a mirror of Chern- 
Simons theory," JHEP 0402, 010 (2004) [hep-th/0211098]. 

[24] N. Drukker and D. Trancanelli, "A Supermatrix model for N=6 super Chern-Simons- 
matter theory," JHEP 1002, 058 (2010) [arXiv:0912.3006 [hep-th]]. 

[25] R. Gopakumar and C. Vafa, "M theory and topological strings-II," hep-th/9812127. 

[26] E. Witten, "Quantum background independence in string theory," hep-th/9306122. 

[27] S. H. Katz, A. Klemm and C. Vafa, "M theory, topological strings and spinning black 
holes," Adv. Theor. Math. Phys. 3, 1445 (1999) [hep-th/9910181]. 

[28] M. Aganagic, M. Marino and C. Vafa, "All loop topological string amplitudes from Chern- 
Simons theory," Commun. Math. Phys. 247, 467 (2004) [hep-th/0206164]. 

[29] S. Bhattacharyya, A. Grassi, M. Marino and A. Sen, "A One-Loop Test of Quantum 
Supergravity," arXiv: 1210.6057 [hep-th]. 

[30] G. Lockhart and C. Vafa, "Superconformal Partition Functions and Non-perturbative 
Topological Strings," arXiv: 1210.5909 [hep-th]. 



37 



Zi(l) = l/4, Zi(2) = 1/(16*-), Zi(3) = (-3 + *)/(64*), Z,(4) = (10 - * 2 )/(1024* 2 ). Z,(5) = (26 + 20* - 9* 2 )/(4096* 2 ). Z,(6) = (78 - 121* 2 + 36* 3 )/(147456* s ), Zi(7) = (-126 + 174* + 193* 2 - 75* 3 )/(196608* 3 ) 
Zi(8) = (876 - 4148* 2 - 2016* 3 + 1053* 4 )/(18874368* 4 ), Zi(9) = (4140 + 8880* - 15348* 2 - 13480* 3 + 5517* 4 )/(75497472* 4 ), Z,(10) = (16860 - 136700* 2 + 190800* 3 + 207413* 4 - 81000*')/(7549747200*') 

Zi(ll) = (-122580 + 381900* + 837300* 2 - 1289300* 3 - 1091439* 4 + 447525*')/(30198988800* 5 ), Zi(12) = (626760 - 8856300* 2 - 18446400ir 3 + 35287138ir" + 30204000*' - 12504375*')/(4348654387200*") 

Zi(13) = (1563480 + 6714000ir - 17252100ir 2 - 40746000ir 3 + 49141894* 4 + 45780780ir 5 - 18083925*")/(5798205849600*") 

Zi(14) = (21382200 - 421152060* 2 + 1918350000* 3 + 2614227910* 4 - 5654854800ir 5 - 3965159223ir s + 1732468500* T )/(3409345039564800* 7 ) 

Zi(15) = (-222059880 + 1271579400* + 3613033620* 2 - 12266517900* 3 - 17757814914* 4 + 28941378130ir 5 + 21727092861ir s - 9162734175* 7 )/(13637380158259200* 7 ) 
Zi(16) = (288454320 - 8196414240* 2 - 54540622080* 3 + 83379537976* 4 + 337956998400*' - 310977507352*' - 354450849984* 7 + 132764935275* 8 )/(872792330128588800* 8 ) 

Zi(17) = (3171011760 + 23555952000* - 71723746080* 2 - 333199608000* 3 + 542885550648* 4 + 1355261623520*' - 1384280129304*' - 1337978574000* 7 + 518021476875* 8 )/(3491169320514355200* 8 ) 

Zi(18) = (4970745360 - 180631896480* 2 + 2270514395520* 3 + 2444801550408* 4 - 18251132155200* 5 - 13590443330584*" + 35949047139936* 7 + 20671882502409* 8 - 9607077219600*')/(377046286615550361600*») 

Zi(19) = (-2636096400 + 24895105200* + 79219113120* 2 - 487774106400* 3 - 852843285000* 4 + 3053792290360*' + 3630439618136*" - 6122444513560* 7 - 4288974330849*" + 1840384320075*')/(55858709128229683200*") 

Zi(20) = (361238119200 - 17220213738000* 2 - 311199827904000* 3 + 324620205099120* 4 + 4438412804448000*' - 2865149087601400*' - 16634573387107200* 7 + 10309541088246858* 8 + 14094700187198400*' 

-4979889726328125* 10 )/(603274058584880578560000» 10 ) 
Zi(21) = (5559930784800 + 65609800200000* - 214704822906000* 2 - 1706798293200000* 3 + 3122732798020080* 4 + 15272968981284000*' - 19948322354095800*' - 51815653479308000* 7 + 47747674792873242* 8 

+ 47396132831956500*" - 18101156239333125* IO )/(2413096234339522314240000* 10 ) 
Zi(22) = (25844898895200 - 1500424524522000* 2 + 46684893872232000* 3 + 35164010074651920* 4 - 740786057283888000* 5 - 402614116951888600* 6 + 3793710085975101600* 7 + 2103302735093392558* 8 - 6256890344362442400*' 

- 3210299116242733125* 10 + 1551114253328062500* n )/(1167938577420328800092160000* 11 ) 

Zi(23) = (-496645968218400 + 7222685259636000* + 23988125937006000* 2 - 241443239284170000* 3 - 452965450178712240* 4 + 2930864234437062000* 5 + 4031113133564569800*' - 14930904869365243000* 7 - 15786905528345919066* 8 

+ 27222797625410117610* 9 + 18233204850331964775*'° - 7915669863827158125* n )/(4671754309681315200368640000* u ) 
Z,(24) = (1204104891528000 - 86634289079553600* 2 - 3819319790489280000* 3 + 2602976790752830800* 4 + 100433912061555072000* 5 - 40665871963026060960*' - 833015249980407648000* 7 + 330789435742955662220* 8 

+ 2433643974077824569600*' - 1163340151816826557764*'° - 1809616125406822207200* 11 + 614736070878601546875* 12 )/(1345465241188218777706168320000* 12 ) 
Z[(25) = (2787296269665600 + 49358959705440000* - 164816144642318400* 2 - 2065776940458000000* 3 + 3923291222794334160* 4 + 32721154659114768000*' - 46361197897271310240*' - 236251940186723032000* 7 

+ 268046321968484065644* 8 + 718649496677642926000* 8 - 615263237145329762916*" - 623636559509281432200*" + 235728636221685481875* 12 )/(597984551639208345647185920000* 12 ) 
Zi(26) = (112416665053492800 - 9539770920889809600* 2 + 679050246669809414400* 3 + 342466609375043833680* 4 - 17962043874947027856000* 5 - 6538704303611999307360** + 175693381923498780781440* 7 + 67906500509898555099452* 8 

- 732942357598839937348800*" - 340632835309262972353644* 1 ' + 1094527855640549701449936*" + 521717090330960815763475* 12 - 258856873289859609765000* 13 )/(3638138012172943574917479137280000 * 13 ) 
Zi(27) = (-2837783495662315200 + 60576452126406216000* + 201778747825504449600*" - 3107566895728375080000* 3 - 5911562189814162401520* 4 + 62328030245503696674000* 5 + 89382763214933595881760*" 

- 596840165742039299652000* 7 - 710285728070533874332788* 8 + 2688501047332900049538860*" + 2631613640469642289362564* 10 - 4607029942628752806862500*" - 2990942157677468764775625* 12 
+ 1309268426599376039671875 * 13 )/(14552552048691774299669916549120000* 13 ) 

Z[(28) = (12969482004169891200 - 1313678612062041825600* 2 - 130700139399240559948800* 3 + 57505552398883116047520* 4 + 5584407193839421531584000*' - 1387144249864003881410160*' - 82911547437043021895278080* 7 
+ 19393348493721346914697928* 8 + 522086937010309473445843200*" - 149459897799711476890866044* 10 - 1306120034852713622360250432*" + 517277319128577543305082654* 12 + 887470472209191527843877600* 13 

- 292941866060855176428388125 * 14 )/(11409200806174351050941214574510080000* 14 ) 

Z[(29) = (356597653882854384000 + 9094131200774299680000* - 30038899805342890862400* 2 - 561876995802308266080000* 3 + 1062648884610900631044000* 4 + 13932441993654741781224000*' - 20002293314753598671462640*' 

- 174087341108384502887088000* 7 + 208625406705598959725168520* 8 + 1102194201726340597528654000*' - 1127922427503232413314046876* 10 - 3106694796352524493370982000*" + 2512129123834667330105351214* 12 
+ 2590737342046361982559266300* 13 - 971173742724208860654943125 *")/(45636803224697404203764858298040320000* 14 ) 

Zi(30) = (1516562211091623696000 - 177090468308691897960000* 2 + 27136431922196521884960000* 3 + 9017847591500109108146400* 4 - 1077113128060105913256480000* 5 - 256681199098892348426646000*' + 17131486119850520142660936000* 7 
+ 4335785790236401386882916920* 8 - 133790555452597949153456208000*" - 42194568890414867791467229900*" + 497101282096648260408501524400*" + 205605092600586074416500476442* 12 - 697543537772745100897507429200* 13 

- 315952706311600959133249618125* 14 + 159786318876990035803260187500 *")/(41073122902227663783388372468236288000000* ls ) 

Zi(31) = (-5491025370413030448000 + 166004632983882092400000* + 541029306660436634040000* 2 - 12171024630433182313800000* 3 - 22742123903644972954327200* 4 + 365896598220831803825700000*' + 521165244159524214764514000*' 

- 5649100224613410554739750000* 7 - 6885817029748362284443224360* 8 + 46819207040545598624314349000*' + 50706666062879865301800728100*" - 194513286041291875665429285500*" - 180463482287010623903232340206* 12 
+ 319234142124859370913055662750* 13 + 202614683135946238835538564375* 14 - 89233886735707696636857328125 *")/(18254721289878961681505943319216128000000* 15 ) 

Zi(32) = (26378192029689051552000 - 3580049562649923242880000* 2 - 757132674251549925288960000* 3 + 214985863491932279351404800* 4 + 48450137344329251746575360000*' - 7380031300608320195089248000*" 

- 1141207573627181061863237376000* 7 + 155907850426836134859071679840* 8 + 12598650074482753279364061696000** - 2023063868341111494879077115200*" - 66659187395884098831853323686400*" 
+ 14997054793962534623552524193424* 12 + 149677247547917349534167645606400* 13 - 51289784809585816523329448010000* 14 - 95114261347756370696738538000000* 15 
+30707052223734699797941597265625* 1 ")/(21029438925940563857094846703736979456000000* 1 ") 

Z!(33) = (938391146991441535392000 + 33397176671939336256000000* - 106969904782051711720320000* 2 - 2869098805616518579104000000* 3 + 5271786469822631334556780800* 4 + 102890141125717999328795520000*' 

- 144488682885133549916265312000*' - 1981982524489386054179544000000* 7 + 2355690689320882868185315724640* 8 + 21374876708906309003768152800000*' - 22562204399707743197488161172800*" 

- 123272970849585219046190213680000*" + 116183975669541868107557485762704* 12 + 327775305305306832840663235099200* 13 - 252942132305492678291871933104400*" - 264752932078174137625582781940000*" 
+98556910328894483095041052265625 *")/(84117755703762255428379386814947917824000000*") 



Z,(34) . (37 18 ll>50:158 10087038880000 - 576507161916513586991232000ir 2 + 181107877834l)55G5974485700( jr" I :194959075;18!)G2450446622752000» 4 - 10032226168135957514;i54545920()[)0>r r ' 4562219208018840347542:1245888(101' 

+ 234599266467263376168123553920000it 7 + 38603239231337392309701577173600^ - 2932684274647031133516687271872000JT 9 - 600271248924227540020652100904640ir 10 + 201599556774660225455882532724800007T 11 

+ 5573654587796424883190789264494160;r 12 - 69709092298931690593778835499708800ir 13 - 26569384887189154463277158157743664ir 14 + 938326778238255483739240232270064009r ls + 40956483693699951777807871375445625;r 16 

- 21010234864640717487636413068500000 ir 17 )/(97240125593549167275206571158[)797930(>4544l)000OO;r 17 ) 

Zi(35) = (-155394189350732306598624000 + 647040692134313988427920000057 + 2029882077930927425094489600057 2 - 6442377169511302()05999072()()000^ 3 - 1159426487470456162673933817600jr 4 + 27184334880173756913462875040000ir s 
+ 374417125018579405893556730304007r fi - 615920597406312013511953377120000tt 7 - 7373945765034185818589543015380807T 8 + 8083939286151916595206710203292000ir' + 889174539914391886768332225865632(V 

- 60943223839020959285920727787464000ir" - 61938023356886006851232342745704688i 12 + 239331994588192630847481782926880400i 13 + 213684392930403426579039783004190352ir" - 380703883988580770014435383951490800»" ^-"^ 

- 2375797756047133726917494201271768751" + 105110725506673319593388326238765625 n J7 )/(388960502374196669100826284632319172018176000000»' 7 ) T5 
Z,(36) - (205449951732963783797952000 - 36007191928401346014826464000ir 2 - 15429603717166905589305894912000»" + 2840446363912624200459178252800i 4 + 1397363237706856558307600471040000ir* - 131357443180004913764646329145600»" 2j 

- 48296798398044665685367034661580800ir 7 + 3883644257708774705265766583576640i 8 + 831449949125205134950635749923584000i* - 74966355578150879844550753874386080ir"' - 7594849087488363684192700143463050240»" 

+ 921494334533293208733930116271590944ir 12 + 35556136643862328691978571274984281600ir" - 6631878927263733277152930079345687248ir" - 73649263521729777257048019164426375424jr" + 22472441260989832462435020189734075550ir" "_p 
+ 44419572428849001959102055192929040000ir 17 - 14087391698607910119184937538270796875 ir' ! )/(168030937025652961051556954961161882311852032000000i' ! ) Eh 
Zi(37) = (344658385327622963268672000+ 166962754916064875607840000007T - 51162197245362439004504736000ir 2 - 1908396918432031940168256000000ir 3 + 3353443144313022700702847884800tt 4 + 93646325534078783973438771840000;r 3 ^ 

- 1260135276885312824332334142144007r fi - 2584302296129015527191445038720000tt 7 + 2947512886694682727712026276981440ir 8 + 42496143753135725779909660878000000ir 9 - 43492244872629368821077100733533920tt 10 "" — ^ 

- 411952846480704125142442G(>9087()720[)Ojr u + 391369473271880097081811721638702944x 12 + 220891258718329683209772117373409440(V 3 - 1942475465461536631413106579427584752;r 14 - 5601595538289187295174864184144292800ir ls ^ 
+ 4154644618393775505851430554713346250ir 16 + 4405698003009276640285707242239312500jt 17 - 1630215031442982020424451720597265625 jr 18 )/(2489347215194858682245288221646842700916326400000(l7r 18 ) ^ 

Zj(38) = (34670979645901746307288128000 - 6793886204041650800608345824000it 2 + 4202171522328588667044126804864000jr 3 + 6023635137638756928290110061952007T 4 - 307858874840519708881339429224960000jt b ! 

- 31540825837298445178949118762681600ir I + 9864416085598242470628358653522201600ir 7 + 1067017024880265019830906086445725760i 8 - 178863182236055742780317653014420096000i» . 

- 23936003744539769855200093232751864480it'° + 1950971531445274034868089224826497687680i" + 350412682409633888332406915496775780576ir 12 - 12409259364678554294486401637878253126400ir' 3 ^ 

- 3139232447115718702513052765143276159248ir 14 + 40954558592324625344372057487446221111488i 1s + 14707382078185621167538839421215582064938i" - 53558352942981242055276392144650651147200ir 17 Q 

- 22737208019032436560743755455032885721875ir' 8 + 11792146368175750691275947360407520187500 i 1 ")/(970546692260171503033792971855671032233257336832000000ir 1 ») ^ 
Z,(39) - (-1805904761210022513031992000000+ 101264771459781609028217668800000» + 302381839841825904791539997088000i 2 - 13176351246073953964062405669600000ir'-22547043029826270566695173241152000» 1 pL, 

+ 744314423610103051803959582360640000tt 5 + 9753403258782605936199595194597312007r fi - 23199798857345391433844810282804640000tt 7 - 26711620734978466698070836639363700800tt 8 E^ 
+ 439479520445053419416040646131248808000jt 9 + 474622609600170794900723628199437876960jt 10 - 51824430882265527831326065668837849960007T 11 - 5352004748469079784950507079841404737440;r 12 

+ 366096053053046440242671995957195795144005T 13 + 35696346279919871803666484520944563156656ir 14 - 138288504259608600641365483955310789896400ir 15 - 120070444515693036962048824855438711760814ir lfi rC^ 
+ 214933454616669430963676986801739145664350JT 17 + 132413932097779200697842535294520418050625ir 18 - 58786281703205223493628273568619944796875 n- lfl )/(3882186769040686012135171887422684128933029347328000000ir 13 ) 
Z,(40) - (3313407099483014501115308160000 - 728706447258906691069678012800000ir 2 - 608182922633989792238430586905600000ir 8 + 73078733921342337274895668091232000ir 4 + 74789865796932340434322044015513600000»' CSTJ JO 
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- 458791053546S4:15970(17I)8101I)38:17918730I1[)0- 1 ° 125699li:i79521422i:i80449250S7677477286272000ir u + 829022951161)5062905991, 50M 1282()572;i0202l)80;;" I lllll.l?»7<) 10l)2HS.54,S.42:174(,:i0l).5Sr,li777;»736S32000ir la 
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- 9776223936589128743853907658094307:i777700(lir" 4287917:127:11)14(1(1872749219.182:117:1545:1:161:1521)11:; |; ' . 68728 1800:1775:i079l,108U55:i(,l)75«9:i«2:ll,l,122276;r" 1 I 8:i:1857:il6l,77012s54670:i57527274757:i(,01(i(,912O0» 17 ZT. 

- 23118022037577776612867744615810440354417500I 18 - 48213251336088521517706522595147840710500000ir" + 15067070464843057084805405574377241085546875 ir 20 )/(3105749415232548809708137509938147303146423477862400000000 it 2 ") 2j 
Z,(41) - (188552554107834346310449599360000+ 12181876275428018759588788800000000ir- 3537742438981452092619102114240000(V - 1791039072418016434645144312800000000ir'+ 2978072491010546421716012885128032000» 4 ^ 

+ 14543842559827389267136693229817600(1000ir 5 - 146933315003754853108030592402866560000tt 6 - 4341672644396621320587818863381536000000;r 7 + 4655916455401019181633684383830965494400JT 8 + 1015908836011553259058831707985665840000007r s ^> 

- 973289650158149361069470922222763358880007T 10 - 1485822441477880712510802305222954666400000ir n + 1333986362874605156945144268101045806838880ir 12 + 13224665210527084220043346689487904097760000ir 13 ^ 

- 11453230318564880605521757863123786256420800;r' 4 - 66752596672330883653121453078480975195280000ir 15 + 552528704282942243018564880425795754228317167r 1B + 16263158196672474301320503518567795319433000(V 7 O 

- 11652874.57026242317272090921421376419326975007T 18 - 125022202563610959798844577657382305614875000jr 19 + 46015866656808256643855214590685881916796875 7r 2n )/(12422997660930195238832550039752589212585693911449600000000 tt 2 ") i^j 
Z!(42) = (652054007140292642485862021760000 - 158609532299761374501366467971200000ir 2 + 186665933069369587242955731029644800000tt 3 + 17667040439233978395848715632434272000;r 4 - 17341808515612802831731351146454137600000ir s CD 

- 1181703323131690599008517421707914880000ir s + 721507120701706831784587777403347345920000i 7 + 52283990298487438398501982107140915894400ir ! - 17688551061484841941422127936622960605440000i° - 158762587308034361184562279658056348190400()»™ qi. 
+ 277866469275880390352640356421283819649856000ir u + 33172764898064489324223178149704271334493280ir 12 - 2797062373428974139738431202487364917367712000ir 13 - 463731667775739696657987276161848827861454400ir 14 

+ 16963673930573775986898603325483482902501433600I 11 + 4037967583397126356247643211196064937066409076i 1 « - 54301142172300990346205630377022886585223939200ir 17 - 18651852806565966140235166219665124918106872500ir' ! E - ' 
+ 69572684576100005245492562401033725156003750000i 19 + 28912911736598093723356727818160605656338671875ir" - 15123842776704225142080715044329580376171875000 ir 21 )/(2191416787388086440130061827012356737100116405979709440000000() it 21 ) 

Z,(43) - (-4690890099680245571830765595520000 + 347706936435165183500269293648000000ir + 980408078432865164250109198857600000ir 2 - 57381500752645218341404421729520000000i 3 - 92542081171004968477059413691146784000I 4 ^ 
+ 44869952713715587354856854916944920000001 s + 5165487961081754007275196622046087040000ir 5 - 171069646711081801655850927133868784000000ir 7 - 1874510351041738546070056952209506174048001 s + 4381198954546784241621461197678965762960000ir' 

+ 4593784730914452699343606984329562512912000ir 10 - 737594206277889342236190840084482525988000001T 11 - 75685608695625981059554748404625147226147360ir 12 + 810639408697052935993089371122080494389820000ir 13 ^ 

+ 8092264491)2 10922182 17.59716 112 17241 108555 17200;. 11 548;<5975170G:I4(>535496:14115!)3:I614025 ll:19772l)l)(H)i- 1 " 5206(>122.529<!0005t)79i)i| 1 ill: Mil '251 i ::) 1 74- r .247:i2ir l,J I 201440976149059211 19():1195:152772:):15()257000207300jt 17 ^ 

+ 47138251755575 HOC h, 1 12 .in.-, in. ^9(11 94 II 50:|2():U!)!)i 1 :i()7()5265195:ll 17950 45 41 190101U9 1 1 1 100751 151)5751)0:.; 13 18707.565500. 197707791. 2 1 1 1 12 . 5i illli " 515:1414204062571 . r-H 

+ 8354300037621534380360695550081236149835546875ir 21 )/(9 7 396301661692730672447192311660299426671M026576486400000000^ 21 ) P 7 -! 

Z!(44) = (142931890873823264893971843552000000 - 385533784091()4()80;l()5189244348472960000jr 2 - 60614534957567443707225912078650112000000jr 3 + 4790543294359345843403437778178955200000jr 4 + 980873638386060440733328907650016616960000(l7r s 

- 360383888599616080396448567562027573792000ir' - 628666211083068466325488711249908753254400000rr 7 + 181373555070592957835428040M757443032800000i 8 + 21457010641852687290301334499478957149834240000ir" 
-6366202949629714859714082118249963841389M400ir' t -43133211761113127020758934306661484113374M80000ir' 1 + 15745638355751087943361823829820040497121592000ir' 2 + 5249902023807224367631030778442404W 

- 270373302125035780750199426666347000631816216480ir' 4 - 37895915015999999636329927324860042910775239296000ir" + 3084963302019749474899609295576812641891628219800ir'«+ 15059428166201725308133960136729014479817495576320()»" 

- 21278839362544516171943770372568304263208365874756ir' 8 - 278275878517177410683024330346304640100600126187200ir'" + 71147751494386429152328551217019038071603976848750ir 2n + 155358605341892695289654301851159713038128800500000ir 21 

- 47950596725606860210962244245014538457634628515625 i 22 )/(169703316015333413923671987883836905721033014479068699033600000000ir 22 ) 



Z 2 (l) = 1/8, Z 2 (2) = l/(32a 2 ), Z 2 (3) = (10 - a 2 )/(512a 2 ) Z 2 (4) = (24 - 32a 2 + 3a 4 )/(49152a 4 ), Z 2 (5) = (1032 - 904a 2 + 81a 4 )/(1179648a 4 ), Z 2 (6) = (360 - 1440a 2 + 2807a 4 - 270a 8 )/(70778880a 8 ) 
Z 2 (7) = (359280 - 918600a 2 + 589034a 4 - 50625a")/(16986931200a 8 ), Z 2 (8) = (20160 - 161280ir 2 + 1366288a 4 - 1228032a 8 + 110565a 8 )/(507343011840a 8 ) 

Z 2 (9) = (145010880 - 757626240a 2 + 1292565904a 4 - 723655312a 8 + 60824925a 8 )/(426168129945600a 8 ), Z 2 (10) = (302400 - 4032000a 2 + 130880400ir 4 - 291792800a 6 + 176521079a 8 - 15025500a 1II )/(1217623228416000a 10 ) 
Z 2 (ll) = (74456323200 - 664676712000a 2 + 2192481033120a 4 - 3105786871600a 6 + 1607656687398a 8 - 133217004375a 10 )/(18410463213649920000a 10 ) 
Z 2 (12) = (558835200 - 11176704000a 2 + 1296197179200a 4 - 5590427673600a 6 + 7745173132312a 8 - 3722093964000a 16 + 303293615625a 12 )/(432032203413651456000a 12 ) 

Z 2 (13) = (194023988928000 - 2657411755891200a 2 + 14431398548212800a 4 - 38219526536737920a 6 + 48507686982371320a 8 - 23823195742030968a 18 + 1954077051898125a 12 )/(5132542576554179297280000a 12 ) 

Z 2 (14) = (108972864000 - 3051240192000ir 2 + 1166600822587200a 4 - 84543478699680007r s + 21741780886785960a 8 - 23339271580526720a 10 + 9641958314010093a 12 - 759070490928750a 14 )/(18871166645108295598080000a 14 ) 

Z 2 (15) = (787945104032486400 - 15409012442388364800ir 2 + 124996926434092352640a 4 - 531883490421787183680a 6 + 1234759142546311157776a 8 - 1454710320020335668472a 18 + 687260380130920046898tt 12 

- 55929616500716083125a 14 )/(2720165444892490160689643520000a 14 ) 

Z 2 (16) = (86306508288000 - 3222109642752000a 2 + 3749940341731584000a 4 - 41498126318134886400a 8 + 174585692946050305920a 8 - 342711141392257397760ir 10 + 312740486692520191712a 12 - 115997092985213588736ir 14 

+ 8880888816551626875a 16 )/(3826166779628997149101916160000a 16 ) 
Z 2 (17) = (1174520269429258752000 - 31153310208889989120000a 2 + 353604351346282465228800a 4 - 2206170741163938861312000a 6 + 8099218379243719626138240a 8 - 17182132282133938856172800a 10 

+ 19149751631150222464481184ir 12 - 8774440277838741548100000a 14 + 709487455486528859765625 a 16 )/(626726118503229733022893867008000000a 16 ) 
Z 2 (18) = (44505389440512000 - 2136258693144576000a 2 + 7045481116676075673600a 4 - 111270721326427856793600a 6 + 701552429702000576219520a 8 - 2223697665491636216762880a 18 + 3677115584554262596762016a 12 

- 2994137064571764051874752a 14 + 1023827859827918055824283a 16 - 76595501984076690615000 a 18 )/(568231680776261624607423772753920000a 18 ) 

Z 2 (19) = (3275305728023106505728000 - 113378060870190461689344000a 2 + 1717995073916883638151475200a 4 - 14779491212904380459162265600a 6 + 78664265877719398398796834560a 8 - 260899435751714064910595652480a 1 ' 
+ 516998276680469455334417352256a 12 - 551137960521620544983472700128a 14 + 246212931459333034614107699550a 16 - 19799702711969287082272171875 a 18 )/(312982009771564902833769014032859136000000a 18 ) 

Z 2 (20) = (101472287924367360000 - 6088337275462041600000a 2 + 53490193257256743955968000tt 4 - 1147279812673273334415360000a 6 + 10180480198229895204569318400a 8 - 47824540468978617999502848000a 10 
+ 126522983630835272734199506560a 12 - 185451912547633642015127859200a 14 + 138679030761041176293852038232a 16 - 44485762576443496864931210400a 18 
+ 3263739072714937244162109375 a 28 )/(414581834294360481313576384601260032000000a 20 ) 

Z 3 (l) = 1/12, Z 3 (2) = (-3 + a)/(48a), Z,(3) = (9 + 108a - 64\ / 3a)/(5184a), Z 3 (4) = (378 + 180ir - 539a 2 + 256x/3a 2 )/(82944a 2 ), Z,(5) = (702 + 6480a - 768V3a - 5701tt 2 + 2304N/3a 2 )/(2985984a 2 ) 
Z 3 (6) = (-2430 + 1998a + 13923tt 2 + 3840V3a 2 - 6673a 3 )/(11943936a 3 ), Z 3 (7) = (17982 + 1487160a - 134784v / 3tt - 235575a 2 - 435456V3V - 2586816a 3 + 1601728V3a 3 )/(3869835264a 8 ) 
Z 3 (8) = (1128492 + 1185840tt - 16264476?r 2 + 3068928v / 3a 2 - 13831992ir 3 + 35414397a- 4 - 17292800V3a 4 )/(123834728448a 4 ) 

Z 3 (9) = (1472580 + 89719488a - 3068928\/3a - 46127556a 2 + 45536256v / 3a 2 - 378069120ir 3 + 52376064v / 3a 8 + 378253327a 4 - 169072128\/3a 4 )/(4458050224128a 4 ), 
Z 3 (10) = (-145260540 + 208275300a + 3644198100a 2 + 632448000\/3a 2 - 2972675700a 3 - 19775216637a 4 - 6435033600V3a 4 + 10570974175a 5 - 335923200v / 3a 5 )/(445805022412800a 3 ) 
Z 3 (ll) = (299312820 + 88250698800a - 2356128000V3a - 12526577100a 2 - 33312384000\/3a 2 - 1023089097600a 3 + 91604563200v / 3a 3 + 78245930511a 4 + 332129548800\/3a 4 + 1443913778700a 5 
-902368702400\/3a 6 )/(48146942420582400a 5 ) 

Z 3 (12) = (8708546520 + 15182154000a - 355130527500a 2 + 53318476800^3a 2 - 643639834800a 3 + 4434569325126a 4 - 1499789952000v / 3a 4 + 4618381766220a 5 - 15721205760073a 5 - 10986823994175a 8 
+ 5450595219200\/3a 6 )/(770351078729318400a 6 ) 

Z 3 (13) = (25651672920 + 6020447083200a - 76624081920^ - 1788699404700a 2 + 3205211212800\/3a 2 - 106194939321600a 3 + 3814758374400v / 3a 3 + 59135521614966a 4 - 61844039654400\/3a 4 + 326960550693360a 5 

- 49162115356416\/3a 5 - 353685210676475a 8 + 164389189305600\/3a 6 )/(83197916502766387200a 6 ) 

Z 3 (14) = (-5599326148920 + 12192598920600a + 327391390155060a 2 + 47611234944000v / 3a 2 - 331983697632300a 3 - 5928588773974566a 4 - 1240653059251200v / 3a 4 + 4694982635549430a 5 - 482647063449600v / 3a 5 

+ 33032941562527173a 6 + 11297648468171520v / 3a 6 - 18749272956026375a 7 + 1083923053516800V3a 7 )/(16306791634542211891200a 7 ) 
Z 3 (15) = (9762153103080 + 8002729127423520a - 80443646277120^ - 842105862572940a 2 - 3248179182904320v / 3a 2 - 251373308727744000a 8 + 6524189288467200\/3a 8 - 14425113463704846a 4 + 102523904471116800V3a 4 

+ 2042292408062643576a 5 - 181471161819298176v / 3a 5 - 62639788677842547a 8 - 699544133395704576V5a 8 - 2646161456861676600a 7 + 1659255242570878400v'3a 7 )/(1761133496530558884249600a 7 ) 
Za(16) = (1131904919350800 + 2906885902032000a - 91901054842198560a 2 + 12485221294694400v / 3a 2 - 296821079176870080a 8 + 2868185584451494440a 4 - 930153171012710400v / 3a 4 + 6507854398692773280a 5 

- 341106145615872000v / 3a 5 - 32541301215263383128a 8 + 14410703617461688320\/3a 6 - 37430599676027727024a 7 + 2308644602346700800 V3a 7 + 85096547210538648225a 8 
-42660063670951193600v / 3a 8 )/(112712543777955768591974400a 8 ) 

Z 3 (17) = (873064837174320 + 610268956142952960a - 3332148259184640V3a - 107766663838898400a 2 + 317593856147742720\/3a 2 - 25458746188524364800a 8 + 327107190466990080v / 3a 8 + 13328971254570106296a 4 

- 14907554349424128000v / 3a 4 + 299477269154064774528a 5 - 11111084392271173632v5a 5 - 189225618460189692840a 8 + 189805167698060630016^" - 802724981961779683200a 7 + 128239157457940337664v5a 7 
+ 914684272980604229075a 8 - 434650602648511872000V5a s )/(4057651576006407669311078400a 8 ) 

Z 3 (18) = (-4357116781128720 + 13351215166858800a + 462502127490651360a 2 + 62013565910016000v^Sa 2 - 406665151540875360a 3 - 16710670410456818376a 4 - 2425393899408936960v'3a 4 + 10314727603578551160a 3 

- 3951295513985433600\/3a 5 + 278995588617672585288a" + 56862983172665456640^" - 213566979579646867848a 7 + 53942665950893260800v / 3a 7 - 1634570622021407629545a 8 - 569046977693061731328v / 3a 8 
+ 965826228941841503075a 8 - 76073671439016652800v / 3a s )/(16230606304025630677244313600a 8 ) 



Z t (L) = 1/16 

Z 4 (2) = (-8 + 7r 2 )/(5127r 2 ) 

Zi(3) = (-8 - 32tt + ll7r 2 )/(81927r 2 ) 

Z 4 (4) = (192 - 560tt 2 - 384tt 3 + 1777r 4 )/(15728647r 4 ) 

Z 4 (5) = (576 + 16896tt - 3120tt 2 - 133767r ;i + 40237r 4 )/(754974727r 4 ) <*6 
Z 4 (6) = (-23040 + 192960tt 2 + 506880tt 3 - 494312tt 4 - 574080r' + 2145157r 6 )/(362387865607r 6 ) || 
Zi(7) = (-115200 - 14653440tt + 1396800tt 2 + 44582400tt 3 - 5718760tt 4 - 252621127r 5 + 72164257r 6 )/(28991029248007r 6 ), 

Z 4 (8) = (6451200 - 106444800tt 2 - 8205926407r :i + 655571840tt 4 + 32062464007T 5 - 16094728007T 6 - 2421603072tt 7 + 8265489757r 8 )/(25975962206208007r 8 ) 
Zi(9) = (135475200 + 57823395840?r - 2912716800tt 2 - 403938877440tt 3 + 24113927040tt 4 + 807646292992tt 5 - 87822722400tt 6 - 394780516992tt 7 
+ 1095693054757r 8 )/(8727923301285888007r 8 ) 



+ 18517120302453483110407T 7 - 237088654639272916480tt 8 - 7449382173461106278400tt 9 + 11631422183825876432007T 10 + 128365438337708670658567T 1 
- 29582301566521923238807T 12 - 6645748001721802646400tt 13 + 2013179571256608129375tt 14 )/(36277729466111972969262612480000 tt 14 ) 
Z 4 (15) = (-574456119164928000 - 4252705107415872307200tt + 37698682820198400000tt 2 + 128051242036867576627200r' - 1090001563845562368000tt 4 



O 



Z t (W) = (-1083801600 + 29578752000tt 2 + 578233958400?r 3 - 328147276800tt 4 - 4901011046400^ + 1776733212800tt 6 + 12068917073920tt 7 - 4240588979256tt 8 

- 7552952651520tt 9 + 24425221038757r 1 ")/(1396467728205742080007r 1 ") 
Z 4 (ll) = (-29262643200 - 35762356224000tt + 981517824000tt 2 + 4561631 13369600tt 3 - 13560288537600tt 4 - 1989404766720000tt 5 + 94612817577600/' + 3271732109716480tt 7 CSJ 

- 321584570195112tt 8 - 1466895830516640tt 9 + 3993861731658757r 10 )/(603274058584880578560007r 10 ) 

Z 4 (12) = (858370867200 - 34978612838400tt 2 - 1573543673856000tt 3 + 607388292403200tt 4 + 23505886701158400tt 5 - 5634782902425600tt 6 - 120944965839667200tt 7 g , — { 

+ 27397474939488192tt 8 + 238604318535249920tt 9 - 659365984391919847T 10 - 1333916246093587207T 11 + 415573168025688757r 12 )/(424704937243755927306240007r 12 ) 
Z 4 (13) = (47210397696000 + 147960574613913600tt - 2277901688832000tt 2 - 30245536120504320007r :i + 47141631161856000tt 4 + 23265188986777436160tt 5 > 

- 5293614954862080007T 6 - 80848090651157913600tt 7 + 3288078092186986560tt 8 + 1182760111459531658247T 9 - 108113785756710703207T 10 - 500486481981871056007T 11 ~Q 

+ 134351991210855206257i- 12 )/(37374034477450521602949120000tt 12 ) ^ 

CD 

Z t (U) = (-1636627120128000 + 93083167457280000tt 2 + 89762748599107584007r :i - 2324234978537472000tt 4 - 209741539423223808000tt 5 + 326199919675069440007T 6 ^ 



- 15394561332688619156275207T 5 + 17858925579392561664000tt 6 + 9220725502522 143368478720?r 7 - 164677338491547053652480tt 8 - 27893283955626979109828608tt 9 & 



3 

bJO 



+ 9889861170124333301904007T 10 + 376907381624044991156935687T 11 - 3251910491479391104317240tt 12 - 152923121003776328285976007T 13 
+ 4059760856437163341839375 7r 14 )/(2037357286816848401953788316876800007r") 
Z 4 (16) = (4595648953319424000 - 347737437467836416000?r 2 - 68043281718653956915200tt 3 + 117911650279097106432007T 4 + 2300873670656175321907200r' 

- 2321025775290727464960007T 6 - 311400506398612521694003207T 7 + 1606592053498785797283840tt 8 + 2110173564767581250204467207T 9 - 122155928015705289203097607T 10 

- 7289117754349459365714657287T 11 + 80055727755935442206522112tt 12 + 11445660423969572006963527687T 13 - 2276118438677390332112630407T 14 

- 5600558998218723099797376007T 15 + 165922990425511426822089375 7r 16 )/(521563465425113190900169809120460800007r 16 ) 



Z 6 (13) = (62293477299179336179200+ 191373715777382807961600V5tt + 53986778775880541529523200tt 2 - 9714396039974354792448000^ 



Z 6 (U) = (75233046700843903872000 + 62316129128911639177996800tt 2 - 179741171859553222100582400V57r 3 - 2920908275659757270602915200tt 4 

+ 7337476756897728787906560000\/3tt 5 + 55771703247774513136798140480tt 6 - 111172680157979357109753200640V57t 7 - 497846711834332956259427062800tt 8 

+ 803056003567155007256530944000V3tt 9 + 19712383739462014769672666482327T 10 - 2654911265312663505530589182976^" - 23497181000006601903858299010787r 1; 

+ 2874604570520265111460706611200\/37r 13 - 1222676176528703161321453271875 7r 14 )/(284929763397262594971958641752329420800007r 14 ) 



Z 6 (l) = 1/24 

Z 6 (2) = (54 - 57r 2 )/(51847r 2 ) 
Z 6 (3) = (189 + 192V37T - 1257r 2 )/(1866247r 2 ) 
Z 6 (4) = (5832 - 8856tt 2 - 13824v/3tt 3 + 84597r 4 )/(1074954247r 4 ) 
Z 6 (5) = (670680 + 1078272^ + 1035504tt 2 - 3623424 v/Itt 3 + 182566l7r 4 )/(696570347527r 4 ) || 
Z 6 (6) = (1574640 + 4461480tt 2 - 31104000\/3V ! - 10628658tt 4 + 71147520\/3tt 5 - 364587057r 6 )/(83588441702407r 6 ) ^ 
Z 6 (7) = (2211581880 + 4551759360\/3tt + 50163219000tt 2 - 45568396800\/3tt 3 - 90128238381tt 4 + 111021346368\/37r 5 - 500749045257r 6 )/(406239826673664007r 6 ) O 
Z 6 (8) = (17856417600 + 539043361920tt 2 - 1638499000320V3tt 3 - 5410591285680tt 4 + 12087329587200v / 3tt 5 + 8110453018896tt 6 - 18418078780416\/3tt 7 a 

+ 87217006850757r 8 )/(363990884699602944007r 8 ) ^ 
Z 6 (9) = (106825423078080 + 260289957150720^3tt + 10951867524084480tt 2 - 5226616796405760V3tt 3 - 71942666862535056+ 1 + 34782176721954816\/37r r> ^ 

+ 819656198967954247T 6 - 73627337521026048\/37r 7 + 311016397666105257r 8 )/(4953187958992196861952007r 8 ) 
Z 6 (10) = (20249177558400 + 2379575970072000tt 2 - 6463906773811200\/3tt 3 - 46213099086722400tt 4 + 99717782505062400\/3+ 5 + 264787782487023600tt 6 

- 459434908272936960V3vr 7 - 356403859223598126tt 8 + 585393319191736320\^7r 9 - 2641995393244663757T 10 ) /(19812751835968787447808000tt 10 ) ° 

Z 6 (H) = (565659007308561600 + 1566936357235200000\/3jt + 184947764670274536000tt 2 - 52720374262219776000\/3tt 3 - 2526481305792155280240+ 1 

+ 667345559627553408000\/37r r> + 99584435555025539868007T 6 - 3679821240531691852800V3tt 7 - 9435237128848062914391tt 8 + 7176870243498043838400VStt 9 13 

> 

- 28986387080182624531257r 10 )/(8666097653052747629671219200007r 10 ) ^ 
Z 6 (12) = (1082521032272064000 + 365904737966443737600tt 2 - 1005513034318442496000\/3r ! - 11574416716220175048000tt 4 + 26694180539585298432000\/3r' § 

+ 1331236833310322383737607T 6 - 245089914171050924236800V37T 7 - 608208196884211255269960tt 8 + 912176891685213084057600V5tt 9 0> 
+ 764080450644042337595544tt 10 - 1051590384346712928115200V37T 11 + 458877788639987997809375tt 12 ) /(610093274774913433128853831680000 it 12 ) J* 



CD 



- 12574989334157367490709212807T 4 + 198742049521922643898613760v / 3tt 5 + 100955036901251531713462387207T 6 - 2018130164120379746906726400V5tt 7 CO 

H 

- 31512503968726555514761885512tt 8 + 9944682708691937722724803584V5tt 9 + 26854902567483870461075407728tt 10 - 18369730532357431360325337600\/3V n ^ 
+ 7175716737492220019983601875 tt 12 ) /(39138703763360246562082231009935360000tt 12 ) £ 



